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Abstract of the Dissertation 

Applications of instantons to hadronic 

processes 

by 

Valeriu loan Zetocha 
Doctor of Philosophy 
in 

Physics 
Stony Brook University 
2004 



Instantons constitute an important part of QCD as they provide 
a way to reach behind the perturbative region. In the introduc- 
tory chapters we present, in the framework of a simple standard 
integral, the ideas that constitute the backbone of instanton com- 
putation. We explain why instantons are crucial for capturing non- 
perturbative aspects of any theory and get a feel for zero mode dif- 
ficulties and moduli space. Within the same setting we explore the 
configuration space further by showing how constrained instantons 
and instanton valleys come into play. 
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We then turn our attention to QCD instantons and briefly show 
the steps to compute the effective lagrangian. We also show how 
single instanton approximation arises and how one can use it to 
evaluate correlation functions. By this we set the stage for the 
main parts of the thesis: computation of r] c decay and evaluation 
of nucleon vector and axial vector couplings. 

Having understood the effective lagrangian we use it as a main 
tool for studying instanton contributions to hadronic decays of 
the scalar glueball, the pseudoscalar charmonium state r] c , and 
the scalar charmonium state Xc- Hadronic decays of the rj c are 
of particular interest. The three main decay channels are KKn, 
r/Ti-K and tj'tttt, each with an unusually large branching ratio ~ 5%. 
On the quark level, all three decays correspond to an instanton 
type vertex (cc)(ss)(dd)(uu). We show that the total decay rate 
into three pseudoscalar mesons can be reproduced using an instan- 
ton size distribution consistent with phenomenology and lattice 
results. Instantons correctly reproduce the ratio B(7nrr])/B(7nrrj') 
but over-predict the ratio B(KKtt)/ B(7nrr](r)')). We consider the 
role of scalar resonances and suggest that the decay mechanism can 
be studied by measuring the angular distribution of decay prod- 
ucts. 

In the next part, motivated by measurements of the flavor singlet 
axial coupling constant of the nucleon in polarized deep inelastic 
scattering we study the contribution of instantons to OZI violation 

iv 



in the axial- vector channel. We consider, in particular, the f\ — a\ 
meson splitting, the flavor singlet and triplet axial coupling of a 
constituent quark, and the axial coupling constant of the nucleon. 
We show that instantons provide a short distance contribution to 
OZI violating correlation functions which is repulsive in the fi me- 
son channel and adds to the flavor singlet three-point function of 
a constituent quark. We also show that the sign of this contribu- 
tion is determined by general arguments similar to the Weingarten 
inequalities. We compute long distance contributions using numer- 
ical simulations of the instanton liquid. We find that the iso-vector 
axial coupling constant of a constituent quark is {g^q = 0.9 and 
that of a nucleon is g\ = 1.28, in good agreement with experiment. 
The flavor singlet coupling of quark is close to one, while that of a 
nucleon is suppressed g\ = 0.8. This number is still significantly 
larger than the experimental value g\ = (0.28 — 0.41). 

Throughout the instanton computation in QCD one employs in- 
tegration over the SU(2) group as the SU(2) parameters are part 
of the moduli space. The techniques of group integration are es- 
pecially useful for computing the effective lagrangian. We there- 
fore present an algorithm for computation of integrals over com- 
pact groups that is both simple and easy to implement. The main 
idea was mentioned before by Michael Creutz but, to our knowl- 
edge, never carried out completely. We exemplify it on integrals 
over SU(N) of type J du(uu^) n , with n = 1,2,3 as well as in- 
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tegrals of adjoint representation matrices J du(R ab ) 
R ab = ±Tr(\ h u\ a u^). 
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Chapter 1 
Introduction 



The physics of fundamental interactions is dominated today by the Standard 
Model(SM), the combined theory of strong, weak and electromagnetic force. 
Since its development in 1970's it has proved extremely powerful in explaining 
the experimental results. The model has been so successful that, for thirty 
years, physicists have been desperately looking for a discrepancy between the 
model and experiment that would point them towards new physics "beyond 
the standard model". 

Quantum chromodynamics(QCD), as a part of the Standard Model de- 
scribes the gluon-mediated strong interactions between quarks. The origins of 
QCD date back to early 1960's. The myriad of observed particles and their 
mass spectrum played then the same role as Mendeleev's periodic table of el- 
ements a century ago: it pointed to the existence of underlying constituents, 
particles that represent a more elementary form of matter. 

In 1964 Gell Mann and Zweig introduced spin-| particles: up, down and 
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strange quarks 1 with fractional charge of | for u quark and — | for d and 
s quarks. Assigning mesons to qq states and baryons to qqq states using 
SU(3) fi avor symmetry led to a good match of known particles and valuable 
predictions of new ones. (Valuable indeed, as Gell Mann was rewarded with 
Nobel Prize in 1969 for his "Eightfold Way"). 

However, the straight-forward quark model had to overcome the difficulty 
of reconciling the Pauli principle and the seemingly quark-exchange symmetric 
function of baryons made of 3 same-flavor quarks with the same spin, like 
A ++ ~ u*u*u' . The solution came with an extra quantum number: color, 
that made it possible to antisymmetrize the baryonic state. Quarks would 
then come in three colors: red, green and blue and the hadrons would be 
composed of white combinations of quarks. Color- ant icolor pairs of quark - 
antiquark would form mesons while baryons would be of type ^q^qflk- 

At this stage, the colored quarks were little more than mathematical ob- 
jects that explained the spectrum of observed hadrons. The unsuccessful 
search for free spin-| particles with fractional charge presented a big obstacle 
for the quark model to become more than a nice mathematical description 
of hadronic spectra. In late 1960's, SLAC-MIT deep inelastic scattering ex- 
periments discovered point-like constituents, "partons", later identified with 
quarks. Eventually the discovery of J/^> in 1974 in both hadroproduction and 
e + e~ annihilation finalized the conclusion that the quarks are real particles 
but are confined to colorless combinations. 

The main conclusion of DIS experiments was the asymptotic freedom of 

proton constituents, which basically means that the interaction of partons is 
charm, bottom and top were added later 
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small for high energy transfer (> 1 GeV). This was an extra feature that any 
theory of strong interactions would need to possess. 

QCD, the gauge theory of quarks based on SU(3) group then emerged 
as the only viable candidate that would incorporate N c = 3 quarks, have a 
different representation for antiquarks and would feature asymptotic freedom 
and confinement. 

On its way to maturity, QCD underwent a long series of experimental 
tests, like ir° — > 77 decay, deep inelastic scattering, e + e~ — > hadrons and 
so on. Most of the early successes were predominantly in the perturbative 
area of large momentum transfer and hence small coupling. In this sector, the 
methods of computation had already been known from QED and the Feynman 
diagrams techniques were readily available. 

A completely different view was provided by lattice gauge theory methods, 
which use discretization of large but finite volume of spacetime and evaluate 
the path integral using Monte Carlo techniques. With increasing computa- 
tional power the lattice QCD is becoming an accessible laboratory for non- 
perturbative physics. The main drawback of lattice QCD is the absence of 
analytic results from which one could get a better insight into the physical 
picture. 

In the non-perturbative regime different analytic approaches have been 
developed that describe the physics at different energy scales. For small mo- 
menta well below 1 GeV the energetically accessible degrees of freedom are 
pions and other low-lying mesons. Therefore effective theories like the ones 
based on chiral lagrangians were natural tools successfully applied to physics 
of low energy pions. On the other side of particle spectrum, one can use the 
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fact, that low energy interactions are less likely to create new heavy quark- 
antiquark pairs due to energy gaps. Therefore heavy quarkonia are accessible 
through non-relativistic quantum mechanical models similar to positronium 
treatment. 

Instantons bridge the gap between the effective field theories and perturba- 
tive QCD. With their origins in the heart of QCD theory, they constitute one 
of the best understood non-perturbative tools. Instantons saturate the U(1)a 
anomaly and provide dynamical chiral symmetry braking. However, they do 
not solve the problem of confinement, which even after 30 years of work is still 
an open question. 

The main tool for instanton computations is the instanton liquid model 
(ILM) developed in 1980's. Is is based on the assumption that the vacuum is 
dominated by instantons. The parameters of the model, the mean density of 
instantons n = 1 fm -4 and the average size of instantons p = | fm were fitted 
from the phenomenological values of quark and gluon condensates. 

Among the successes of ILM are calculations of hadronic correlation func- 
tions, hadronic masses and coupling constants. A better understanding of 
chiral phase transition has also been achieved by studying instantons in finite 
temperature QCD. 

In addition to that, we would like to identify direct instanton contribu- 
tions to hadronic processes. One possibility is to go to very high energy, for 
example in DIS. In this case, instantons are very rare, but they lead to spe- 
cial processes with multi-gluon and quark emission, analogous to the baryon 
number violating instanton process in the electroweak sector. 

In this work we focus on another possibility. The instanton induced inter- 



1 



action has very peculiar spin and flavor correlations that distinguish instantons 
from perturbative forces. We study two systems in which unusual flavor and 
spin effects have been observed: the decay of r\ c and the so-called "proton spin 
crisis" . 

First part of this work analysis the decay of r\ c . The dominant decay chan- 
nels have a specific structure of final products in which all pairs of constituent 
quarks are present. The resulting vertex is of type (cc)(ss)(dd)(uu) which 
points to a possible instanton-induced mechanism behind the process. 

The experimental evidence that quarks only carry about 30% of the proton 
spin is contrary to the naive quark model predictions and it implies a large 
amount of OZI violation in the flavor singlet axial vector channel. The second 
part of this work studies the instanton contribution to this process. 

1.1 High school instantons 
1.1.1 Toy model instantons 

Instanton is a finite-action solution of the euclidean equation of motion. Since 
it's discovery in QCD in 1975 by Belavin et al [1] it has been enjoying con- 
siderable attention. In this section we will try to convey a feel for instanton 
physics without dipping into the details of any particular theory. 

Many of the features of instanton computation can be easily explained on 
trivial examples that do not require more than high school integration. We 
will use the standard integral as an analog for the path integral. This way the 
ideas of the computation will be unveiled in full light without the technical 
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difficulties blocking the view. 

The single most important object in quantum field theory is the partition 
function, which can be represented as the path integral 2 

Z[J] = J V<S>e~ s[ *' J] 

where is the field and S[$, J] is the action that also depends on the source 
J. Any correlation function can be computed once we know the partition 
function. Schematically: 

< $ x $ 2 >= Sj^ZiJ] 

where J%, J 2 are the sources corresponding to $1, $2- 

It is usually not a trivial exercise to compute Z[J\. In fact, most of the 
time we are forced to rely on some kind of approximation. The standard 
technique is to expand the action around the trivial minimum ($ = 0), keep 
the quadratic terms and treat the rest as a perturbation: 

e Sm = e -S[0]-|*S"[0]*+O(*3) = e -i*S''[0]*^ + + . . . ) ^-g 

However, the perturbation approach usually does not, and can not, reach 

all the 'dark corners' of the theory. The reason is the possible existence of 

other, non-trivial minima of the action, located in a different sector of the 

configuration space. To get a better feel about the approximations to the 

path integral, let us consider the most trivial toy-model for the path integral: 
2 For simplicity we take directly the euclidcan spacetime 
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Figure 1.1: Action of a tilted double well potential. 



QFT for a field in one point in spacetime. The configuration space is now 
considerably shrunk to one single variable, and the path integral is just a 
standard integral. As an example of theory with 'instanton', let us take a 
tilted double well potential shown in Fig. II. H 

= $ 2 (6 2 ($-a) 2 + m). (1.2) 



The partition function is: 

/+oo 
d<$>e~ s W (1.3) 
■oo 

Fig. 11.21 presents the graph of e~ s . The goal is to compute the area under 
the curve. The approach corresponding to standard perturbation theory is to 
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Figure 1.2: Exp(-S) for a tilted double well potential. 



keep the quadratic terms in exponent and expand the rest: 



e 



= e 



(m+a- 



2 6 2)$2 



[l + 0($ 3 ) + ---] 



(1.4) 



The graphs of zeroth and first order expansions are shown in Fig. 11.31 
It is clear that in order to capture the full integral one needs a large order 
expansion. 3 We can get a much better approximation of partition function by 
including the 'bump' from the very beginning: just compute the additional 
area by repeating the perturbation approach for the smaller bump. Mathe- 
matically, its peak is at the minimum of the action, which is nothing else than 
the instanton. This separate term is the 'instanton' contribution: 



3 As we will see later, in QCD the graph under the smaller 'bump' is not accessible at all 
this way, as the bump happens in a disconnected part of configuration space. 




(1.5) 
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Figure 1.3: Zeroth (dashed line) and first (dotted line) order expansion around 
the trivial minima of the action are far away from giving a reasonable approx- 
imation to the exp(-S). Bad news for perturbation theory: with higher order 
expansions, it gets worse before it gets better. 



The total area under the curve would be the sum of the results of perturbation 
approach under both total and local minima of the action: 

Z = Zq + Z^ 

The zeroth order expansions under both minima are shown in Fig. 11.41 It is 
again obvious that this is a much better approximation of the integral than 
what we started with. A choice of a = 1, b = 10 and m = 0.4 gives Z = 0.55, 
Zj = 0.43 and Z + Zj = 0.98. The numerical integration gives Z = 0.99. 

In real models trivial fields and instantons live in separate, disconnected 
parts of configuration space. Therefore one not only needs to add instantons 
to get a better approximation with less computation, but has to account for 
instanton part in order to achieve the right result. 
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Figure 1.4: The combined contribution of zeroth order expansions around 
both minima gives a clearly better approximation to the full integral. 

1.1.2 Zero modes and moduli space 

There is a long shot from the trivial toy model to real world physics and many 
difficulties arise on the way. One of them, the zero modes and the moduli 
space is an omni-present feature that has to be dealt with in every instanton 
computation. The idea is easily explained on a 'next-to-trivial' model of two 
spacetime points, i.e. a two-dimensional integral. 

Let us consider the following 'sombrero' action, depicted in Fig. 11.51 with 
the graph of exp(-S) in Fig. 11.(11 

As a result of the presence of rotational symmetry, we have a continuum of 
minima of action. The set of all 'instantons' is a circle with radius 1. The 
parameters that describe it are called collective coordinates. In our case we 
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Figure 1.5: Sombrero action. 



can use an angle G (0, 2tt): 

{x,y)i((p) = (cos 0, sin 0). 

The parameter space of all instantons is called the moduli space and in our 
case has the topology of a circle. 

Suppose we did not have the err function available and decided to compute 
the integral numerically in a similar way we did it in the previous case of 
one-dimensional integral. We would need to expand around the minimum 
of the action. But which one should we choose? Let's say we just pick an 
arbitrary minimum. The zeroth order expansion will now have a form of a 
two-dimensional distorted Gaussian bell, with different curvatures in different 
directions. If we press ahead and try to compute the volume under such a bell, 
we quickly run into difficulties: the result is infinite. Obviously this is not a 
problem of the integral, but of the method we used: the radius in one direction 
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Figure 1.6: Exp(-S) for the sombrero action. 



of the approximating two-dimensional Gaussian bell is infinity(the surface is 
flat in one direction corresponding to the symmetry direction). Therefore none 
of the minima alone is suited for an expansion. 

This kind of difficulty appears every time we deal with a symmetry of the 
action that is broken by instanton. For each symmetry, there is an associated 
direction in the configuration space along which the action does not change. 
The vector that points in this direction is the zero mode. In our case it is the 
vector corresponding to a rotation, i.e. d^, with <f) being the azimuthal angle 4 : 

= I4 ^' V ^ ^ = Sin ^' C ° S ^ 

The solution of the problem is obvious: one needs to turn the symmetry into 

an advantage, not a liability. The change of variables to polar coordinates is in 

4 In general one obtains the zero mode by differentiating the instanton solution with 
respect to the collective coordinate 
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Figure 1.7: At any point (x*,y*) there is a zero mode direction and a 'per- 
turbative' direction. We integrate perturbatively along the ray while the in- 
tegration over the azimuthal angle just gives the volume 2ir of the moduli 
space. 

order. To use the lingo of instanton physics, we will integrate over the direction 
of the zero mode non-perturbatively, leaving the perturbation method for the 
modes perpendicular to this. In other words, in one direction nothing changes, 
so we should just compute the one dimensional integral in the perpendicular 
direction and then multiply the result by the volume of the space in the zero 
mode direction. 

Let us take as a starting point of our integration an arbitrary point (x*, y*) ^ 
(0, 0). The one-dimensional space perpendicular to the zero mode in this point 
is the line (x, y) ■ = — xsin</> + ycos<j) = with tan0 = - see Fig 11.71 
For every point on the line, we only want to integrate radially, and leave the 
perpendicular direction for later. This can be achieved by introducing the 
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following delta function in the integral: 



S((x,y) ■ 2$) = S(— xsin0 + ycoscf)) 

which forces the integration points lay on the line. To obtain the contribution 
over the whole plane we only need to integrate over the angle, with a weight 
that makes the whole insertion a unity 5 : 



f-2-K 

A(x,y) / 6(— xsincj) + y cos d>)dd> — 1 
Jo 



It is easy to compute A(x,y) = \\J x 2 + y 2 . The whole integral then is: 

J dxdy e~ s ( x,v ^-\/~x 2 + y 2 J 5(— xsin0 + y cos< 
After a change of variables 



x = x cos 4> — y sin i 



y' = —x sin + y cos t 



and evaluation of integral over y' one obtains 



2tt i ;>+oo 

6<f>x- dx' e-^'^Vx 12 (1.6) 



The result is nothing else but the volume of moduli space multiplied by the 

integral in radial direction (integral over non-zero modes). Now we could ex- 
5 This is nothing else than the famous Faddeev-Popov unity insertion 
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pand the radial integral around the 2 bumps it contains in the very same way 
we did it in the example from the beginning of this section. 

The main step in dealing with zero modes is their separation from 'pertur- 
bative' direction by requiring the scalar product to be zero: 

{x,y)-% = (1.7) 

Another interesting interpretation is the following. The question is around 
which of the continuum of instantons one should expand. A very intuitive ap- 
proach would be to expand different parts of space around different instanton, 
mainly, the closest one. Let us take again the starting point (x*,y*). The 
closest instanton is given by the minimization of the distance 

min^lKx* , y*) — (cos0,sin0|| 2 

A differentiation w.r.t. gives the same condition 11.71 One would then sep- 
arate the configuration space into blocks, each one of them being dominated 
by the closest instanton. In our case a block would be a ray with angle <j) 
and the points on the ray would be under the 'jurisdiction' of the instanton 
at (cos 0, sin 0). Adding the contributions of all rays would lead to the same 
result 11.61 

1.1.3 Valley instantons and all that 

We should be on a pretty good footing now that we know how to deal with 
a continuum of minima. For any integral we would find all minima and ex- 
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pand around them, paying special attention to treatment of zero modes. This 
kind of 'turning the crank' could actually lead to a poor result in some cases, 
when there are a lot of important points without satisfying the condition of 
minimum. 

The presence of an 'almost zero mode', a direction with very low gradient, 
would make the Gaussian bell around the minima a poor approximation of the 
integral. To get a better understanding, let us take again the sombrero action 
and slightly tilt it: 

S = (x 2 + y 2 - l) 2 - ex 

The graph of the tilted sombrero action with e = | is shown in Fig. 11.81 The 
surface now exhibits a single global minimum - the instanton - close to the 
point (1,0) and a saddle point close to (—1,0). The exponent of tilted action 
is shown in Fig. 11.91 

If e is small, the volume under the surface of e~ s of the tilted sombrero will 
not differ significantly from the original one. In fact, one could obtain an e- 
expansion of the integral. However, our method of expanding around minima 
would fail, as the single Gaussian bell centered at the global minimum of S 
would not describe well the whole surface. 

Since e is small, we should, in principle do something similar to e = case 
discussed in the previous section. For that reason, we first need to identify 
the important points - 'valley instantons' - that are dominant, in some way, 
for exp(— S) in their restricted vicinity. Let us postpone for a moment the 
strict definition of the valley points. Intuitively it is clear, that in our case, 
the special points would form a loop e-close to the former circle of instantons. 
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Figure 1.8: Tilted sombrero action. 



Once we have the valley trajectory, we would integrate perturbatively in the 
sector perpendicular to valley in the same way we treated the non-zero modes 
before. Now the final integral along the valley trajectory would not simply 
give the volume of the would-be moduli space, as the result of integration 
over non-valley directions depends on the parameter of the given point on 
trajectory. 

We still have not defined precisely the valley trajectory. There are actually 
two approaches with slightly different results: the streamline [3, 4] and the 
so-called proper valley [5] method. 

The streamline method requires a starting point located on the trajectory 
but different from the global minimum. In our case it could be the saddle point. 
Having the starting point, the instanton valley is constructed dynamically by 
following the highest slope downwards. It is exactly the path a stream of water 
would follow - hence the name. 
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Figure 1.9: exp(-S) for tilted sombrero action. 



The proper valley method features more stability and does not require any 
starting point. The valley is now made of points with the lowest gradient 
along the contour of equal action. One should therefore draw the contour 
lines ('izoactas') and for each of them identify the point with lowest gradient. 
Joining these points would give the trajectory. 

Once we have the trajectory, we can start turning the crank again: at each 
point of the valley integrate perturbatively over modes perpendicular to the 
trajectory and then integrate over the trajectory of the valley. 

Another approach to the same problem uses so-called 'constrained instan- 
tons' [6,7]. The idea is very simple: slice the configuration space in surfaces 
given by a family of functions. In our case one could take the slices generated 
by the family of vertical planes x = a, with a G (—00, +00). For example, 
the slice generated by x = is just the standard double well potential in y 
direction. Each slice would feature constrained minima which one could use 
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to compute the integral over the slice. At the end, one would only need to 
sum over the slices, i.e. integrate over a. 

The ambiguity of choosing the slicing function makes this approach less 
appealing. One has to have some physical intuition to use it with success. 
In our case, the almost-symmetry of the graph would point to slicing by rays 
—x sin + y cos = 0. 

There is a lot one can learn about the methods of path integration from 
a standard 2 dimensional high-school integral. We will now take the earned 
intuition and apply it to instantons in QCD. 



In this section we will focus on providing the main ideas for the derivation 
of the effective Lagrangian as well as on setting the stage for using the sin- 
gle instanton approximation for computing correlation functions. A thorough 
review on instantons in QCD can be found in [24]. 

As mentioned before, instantons are finite action solutions of equation of 
motion in Euclidean spacetime. For a Yang-Mills theory, the partition function 
reads: 



and F^v = d^A u — d u A^ + [A^, A u ] is the field strength. Let us for simplicity 
consider SU(2) YM theory. The requirement of a finite action leads to fields 



1.2 Instantons in QCD 



z 




where 
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spatial S group S 




Figure 1.10: The instanton at infinity is a mapping from spatial 3-sphere to 
3-sphere of SU(2) parameters. Each such mapping is characterized by integer 
winding number. 

that tend to pure gauge at infinity: — > u~ 1 d^ 1 u, with u G SU{2). The 
infinity is topologically a 3-sphere. The instanton at infinity is therefore a 
map from the spatial S 3 to S 3 of SU (2) parameters. Every such a mapping 
is characterized by a winding number, an integer that shows how many times 
one sphere is wrapped around the other by the map. 

The whole configuration space of finite-action gauge fields is then separated 
into distinct sectors characterized by different values of winding number. Figs. 
II. Ill and II. 121 show schematically the action S and e~ s in the sectors of n = 
(zero winding number), n = 1 (one instanton) and n = 2 (two instantons). 
The sectors are completely separated by infinite action walls, in the sense that 
any trajectory from one sector to another will have points where Sym — oo. 
It is exactly because of this separation that one has to account for instanton 
contributions, as one can not retrieve them from n = sector. 

Before plunging into any instanton computation, one has to deal first with 
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Figure 1.11: Schematic representation of quadratic approximation in different 
sectors of the configuration space. 



gauge invariance. As any other symmetry, it induces zero modes and flat di- 
rections, derailing the perturbative approach. The way around this difficulty 
was explained on the toy model of section ITTT1 integrate over zero modes di- 
rection non-perturbatively, then compute quantum fluctuations perpendicular 
to zero modes. In the case of gauge invariance, this means fixing the space for 
quantum fields by requiring 

K A T = 

where is the covariant derivative in a background field. After exponentiat- 
ing the Fadeev-Popov determinant, the action becomes a functional of ghost 
fields b and c: 

Z = JvAVbVc e - s[AAc] 



21 



1.2.1 Effective Lagrangian for gauge fields 



The instanton solution of the equation of motion reads: 




2P 2 



— u — u 

(x - x ) 2 ((x - x ) 2 + p 2 ) 2 



V 







(1.8) 



where is the instanton configuration with the center at Xq and orien- 
tation given by the SU(2) matrix u. Here r a is the Pauli matrix and f}^ v is 
the t'Hooft tensor with properties given e.g. in [24]. 

Let us now explore the semiclassical approach, i.e. perturbation theory 
around the minima of the action. Expanding the action in the functional 
Taylor series and keeping the terms up to second order we obtain, for the 
gauge sector: 



Even after fixing the gauge, we are still left with a rigid gauge symmetry 
that brings AN — 5 zero modes, where N is the number of colors. Besides 
that, there are 5 flat directions corresponding to scaling and translational 
invariance. The remaining 5*0(4) rotational symmetry does not bring any 
new zero modes, as any such rotation can be undone by a gauge rotation. In 
other words, the SO (A) rotational zero mode points to the space unavailable to 
quantum fluctuations, as they were made perpendicular to gauge zero modes 
by fixing the gauge. Therefore altogether there are 4iVzero modes. 

The non-perturbative integration over the zero mode directions and subse- 
quent computation of quantum fluctuations limited to Gaussian approximation 



e -S _ e -S[A^}-\AM A A+... 
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Figure 1.12: YM action in different winding number sectors 
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X 



det' M A \ 1/2 ( det'Mgh 



reg 

where prime on determinants means that zero eigenvalues are not included. 
The computation of determinants was performed in [2] and the full result is: 



Z f d f 1 

— = C J d 4 xo-^r J dtt-^e"^ = I dji! 



For completeness, the constant is 



(N- l)!(iV-2)! 



C = - e -«(l)-2(AT-2)a(i) 



with a(l) = 0.443 and a{\) = 0.145 
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Let us now push the computation one step further and obtain the instanton- 
induced effective Lagrangian. We are interested in computing the correlations 
of the type: 

< Al{x)A b v {y) > I= I JvAVbVce- s ^A«(x)Al(y) 
Expanding near instanton A^ = A^ + Aft we obtain: 



< A«(x)A b u (y) >j= df,jAf d \x)A b / d \y)+ < A^ u \x) A b ^ u \y) > 



Instead of computing the correlations this way, we are interested in an effective 
potential V/ such that for large distance, \x — Xq\ 2 » p 2 , one retrieves the 
instanton fields: 

< Al{x)Al{y) >j=< Al{x)A b u {y)Vi > pe n 
One can check, that Callan Dashen and Gross (CDG) potential [26]: 



V/ = V C dg = / dji! exp 



r-1 „2 



9 



provides exactly this. To see that this is indeed the case, let us compute the 
gauge field propagator in the instanton background for large distance. The 
second order expansion in rjF gives, schematically: 

< A M (x)A,(y)V 7 >= J dm < A^x)A u (y)r]FriF > 
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Disregarding the vacuum bubbles and next order corrections in the coupling 
constant, one is left with terms like 



which is just the large distance limit of instanton field. 

The anti-instanton fields lead to similar expressions with the substitutions 
u u , 77 <-> fj . 

1.2.2 Fermions and t'Hooft vertex 

We will introduce now the fermionic degrees of freedom and study how this 
affects the computation in the instanton background. 

t'Hooft's discovery of left-handed zero mode of the Dirac operator in the 
presence of instanton came a little bit as a surprise with huge implications. 
First of all the zero mode renders the tunneling amplitude zero in case of 
massless quarks, as the tunneling is proportional to the determinant of Dirac 
operator. However, the instanton contribution to some Green's functions is 
non-zero and clearly distinguishable from the perturbative contribution. In 
these correlations the small mass parameter m from the Dirac determinant 
cancels against — from the zero mode propagator, as we will show below. 

The specific helicity of zero mode also pointed to chiral symmetry braking, 



< A-(x)F^(x ) > 



(2tt)2 ( X - XoY 



This way every A^ix) field will contribute a factor of 
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and in fact solved the r] puzzle. 

Let us now dip into some details of t'Hooft's computation. The fermionic 
part of QCD action reads: 

S[i',J]=J dxi>(x)(p + m)ip{x) - JJ dxdy tJj s (x)J s t(x,y)i) t (y) 

where J s t is the source used to generate bilinear fermionic fields. 

The massless Dirac operator in the instanton background has a left-handed 
zero mode: 

^( x ) = E. I £l±^$ 

with $ am = ^=-. The semiclassical result will now have an extra determinant: 

£ = J dfilp N fm N fe [-lN f H l M ) )p+2N fa (m] = J dfiIJ 

This shows that the tunneling is suppressed in the presence of light fermions 
by a factor of m f. 

To compute the Green's functions we differentiate with respect to the 
source: 

< VcfaWpiv) >-- 



5 Z[J] 



5J$,(z,y) Z[0] 
Z[J] depends on J through detM^: 



j=o 



det = det(-p5 st + J t 



St) 



Let us take for simplicity the propagator for one flavor. Then the deter- 
minant of Dirac operator with the source is given in terms of the eigenvalues 
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det My, = Y[ X J n , 

n 

where satisfy (—p + J)ip = Xfy- Expanding \„ in powers of J we get: 



where a n ■ J = jjdxdy a n (x, y) J(x, y). The propagator then involves: 



S J af 3 [(«o • J)(Ai +ai-J) 



l j=o 



j=o 



The final evaluation at J = renders all but one term zero. After simplification 
with Z[0] we obtain 

m 

which is nothing else but the first correction to zero energy due to source, i.e. 
standard perturbation theory 



(a -J) =<V°W°> 



Then the zero mode part of the propagator reads: 

< Mx)Mv) >= ^°«0»0i$(i/) • 

Note the — factor that will cancel with m from the Dirac determinant to give 
finite correlation functions. 

In the case of many massless quarks, the operator (—p5 st +J s t) has Nf times 
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degenerate ground state. The surviving term in a Green's function is then the 
product of perturbed eigenvalues, which can be written as the determinant of 
'perturbation matrix': 

~[{a& ■ J) = det st < J st \^° >= det st Idxdy ^°*(x) J st (x, y)ip°(y) . 
i=i J 

This clearly points to the following properties of the Green's functions in the 
zero mode approximation of the instanton computation: 

• at least 2Nf fermions participate (otherwise the contribution is of the 
higher order than zero modes). 

• It has a determinantal structure in flavor, all flavors must be there in 
pairs 

• Quarks propagate in the left-handed mode, while anti-quarks in the 
right-handed mode. 

• One can not have 2 fermions propagating in the same zero mode( this is 
the instantonic version of Pauli principle) 

The graphical representation of the instanton vertex is shown in the Fig. 11.131 
Let us now construct the fermionic effective lagrangian. The idea is the same 
as for the gauge field effective vertex: we need to find L t i H such that: 

< ij a (x)lpp(y)L tl H >pert= $ljj)inst = ~ ifjx)^ fi{y) 



P 2 {V-to^fyt 



1+75 



rmr z 



[x - x 0j 



- fa 



1.9) 
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Figure 1.13: Instanton vertex contains all the pairs of light quarks. Particles 
are left-handed while anti-particles are right handed. 

t'Hooft proposed a Lagrangian of the type: 

S t 'H = Jd[i I jK(tpu)(iJtp) 

with K a constant and u a spinor to be determined. First order expansion of 
the exponent of action gives for fermion propagator: 

< u*)Mv)(MW) >= { 2 J ( ~ t y) G^t'-^L 

which leads to the r.h.s. of II. 91 if we choose tu = — , K = Ap 2 n 2 /m 

To make it gauge symmetric we average over all possible orientations of uj 
to get: 

Combining the many-flavors fermionic and gluonic effective lagrangians, 
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the large-distance effects of instantons can be represented by the following 
lagrangian [26-29]: 

Cj = | J] m q p- 27r 2 p 3 q R (ui 2 tf + l -t a R aa 'fj a ; u a^ q L 

x exp f-^pifj^G^A dz^-dp dU. (1.10) 

where t a = \\ a with tr[A a A fe ] = 25 ab are SU(3) generators, 1 2 = diag(l, 1,0), 
rf^ is the 't Hooft symbol and a^ v = |[7 M , 7^]. The instanton is characterized 
by 4N C collective coordinates, the instanton position z, the instanton size p, 
and the color orientation U G SU (N c ). We also define the rotation matrix R ab 
by R aa 'X a ' = UX a W. For an anti-instanton we have to replace L R and 
fj *-* r]. The semi-classical instanton density d(p) is given by 



N f /o 9\2iV, 



„ N d Q (p) 0.466 exp(-1.679iV c )1.34 JV / /8tt 

dM = — ^ w-i)iw-2)i "~ exp 



2 



57T 



<?(p) 2 J 



(1.11) 

where g(p) is the running coupling constant. For small p we have d(p) ~ p 5 ~ 5 
where b = (lliV c )/3 — (2Nf)/3 is the first coefficient of the beta function. 

The need for this lagrangian is not obvious from computations of correla- 
tion functions, as one can do well without it (better actually), since direct use 
of instanton solution and fermionic zero mode is correct at any distance. How- 
ever, the effective lagrangian becomes a great tool if one tries to compute the 
instanton contribution to matrix elements of type < r] c \KKn >. Our inability 
to write pion fields in terms of more fundamental quark creation/annihilation 
operators renders the straight-forward approach of computing correlation func- 
tions inapplicable. As we shall see in chapter |2l there are still ways to employ 
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the effective lagrangian and compute the matrix element as 

<^KK*> in , t =< lc \J d » IJ V t , H V°°°\KK.>. 

Another great advantage of the effective lagrangian is that it displays trans- 
parently the physical properties of the interaction. This way one easily gains 
physical intuition by thinking in terms of Feynman diagrams. 

1.2.3 Dilute gas approximation 

The instanton is a theoretically tractable non-perturbative object. However, 
the QCD vacuum proved to be much more complicated. A good descrip- 
tion has been achieved by constructing the vacuum from instantons and anti- 
instantons. Within this framework, dilute gas approximation leads to a man- 
ageable computation, where multiple instanton effects can be computed from 
single instanton. In this section, following [72], [73] we will present the main 
ideas of dilute gas and single instanton approximation. 

Configuration space of finite-action gauge fields consists of disjunct sub- 
spaces characterized by different winding numbers. The vacuum expectation 
value of an operator II is given as the sum over these sectors of different ho- 
motopy number n: 

< " 0C|n|MC> = s^;[ DW i„^» (112) 

Semiclassical approximation amounts to getting the minimum of action in 
each sector (solution of the equations of motion) and performing the Gaussian 
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approximation around this solution. In each sector, a multi-instanton config- 
uration is the true minimum of the action. However, it has been argued that 
the superposition of well separated n + instantons and n_ anti-instantons, such 
that n + — n_ = n has a much higher entropy then the true minima for the 
sector with winding number n, therefore dominating the path integral. . The 
expansion around these approximate solutions provide a very good description 
of the true vacuum of QCD. Moreover, the Gaussian approximation proves to 
be expressible through the measure of a single instanton, in case one neglects 
the interaction between the instantons. Neglecting other minima, the path 
integral (jl.l2j) can therefore be written as a dilute gas: 



where d/j,± is instanton measure of moduli space. Phenomenological estimates 
give the density of instantons to be of the order of n ~ 1 fm -4 while the mean 
instanton size has been found to be p ~ 0.3 fm. The dimensionless parameter 
p 4 n — 0.008 shows that the instanton liquid is dilute and therefore one can take 
it as a parameter for expansion of the path integral. The same parameters show 
that for distances <C lfm it is therefore reasonable to expect the one instanton 
to provide the dominant contribution. One should not, however, that small 
distance requirement and diluteness of the liquid are two independent aspects. 
The diluteness renders the expansion in density of instantons meaningful, while 
the small distance justifies the use of a single instanton approximation. 



E 



i 



f(dn + r+(dn-) n -(n) A 
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Neglecting the second and higher orders of the density, one arrives at: 

(vac\U\vac) = (U) + ^ I d/x±[(n)± - (n) ] + 0(ji%) = (n) + SU (1.14) 

+- ^ 

The above expression represents the single instanton approximation (SIA) and 
we will be using it throughout the present work 6 . 

1.2.4 SIA versus t'Hooft effective lagrangian 

Before tackling correlators in the SIA approach, let us comment on some 
generic features of the calculations and highlight the link to the t'Hooft ef- 
fective lagrangian. 

As the name points out, the effective lagrangian induced by instantons is 

2 

only valid for large distances compared to the width of instanton: ^ < 1. 
With the typical size of the order of .3fm, the above condition is satisfied with 
~ 10% accuracy already for distances of ~ lfm. One is therefore endowed with 
an additional way of computing large distance correlators. More importantly, 
the effective lagrangian provides an intuition as to what diagrams are dominant 
and which ones disappear completely. The purpose of this small section is to 
show how well-known characteristics of t'Hooft lagrangian are reflected in SIA 
approach. 

Let us then recall the main features of effective vertex. First of all, it is 

based on zero modes only, so there are always additional terms besides the 

6 Let us briefly comment on SIA versus other first order dilute gas approximation. SIA 
does not account for any correlation between instantons. This can be achieved, simply 
speaking, by using an effective mass m* instead of the current quark mass m, with the 
instanton interactions being funnclcd into the value of to*. For more details see [8] 
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ones given by t'Hooft Lagrangian. However, once the zero modes give a non- 
zero contribution, all the non-zero mode terms are suppressed by the quark 
mass m and can be neglected in the m — > limit. 

One can read all the fundamental features of the zero mode propagation 
from the form of the lagrangian: YludsQR (l + | <jM ^ ? 7^ cra ) Ql- All massless 
quarks have to participate, they come in q^q^ pairs and their helicity flips. The 
determinantal structure also restrains the species of quarks from participating 
with more than one pair (Pauli principle). 

None of the above rules have to be imposed by hand in SIA - they are 
already incorporated by means of chirality of the zero mode or rules of Wick 
contractions. It is instructive to see how it works on some simple examples. 

Helicity of the quarks flips through the instanton vertex. A vector insertion 
at point y into a quark propagator from x to z gives rise to 

S ZM (x,y)YS ZM (y,z) = ~M^Uv)Y ~Mv)^ \z) (1.15) 

m m 

The chirality of the zero mode *$?o(y) = 7±^ / o(y) an d anticommutation relation 
{7 / \7 5 } = renders the diagram zero since 

vt(y)-fMy) = ^o(y)i±i fl i±^o{y) = ^t(y)YiTi±^o{y) = o 

All massless quarks participate - is just a statement that, whenever possi- 
ble, the zero mode propagation is favored to non-zero mode due to — factor. 

The most interesting feature of t'Hooft Lagrangian is the Pauli Principle: 
no 2 identical quarks can propagate in the same state. This stems of course 
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from anticommutation of fermion operators. But that is true independent 
on the species. What makes it work is that whenever there are 2 identical 
quarks propagating in the zero mode, there is an additional diagram that gives 
exactly the same contribution with opposite sign. To illustrate this, consider 
the correlation of scalar operator uu: 



(uu(x)uu(y)) = —Tr[S(x, y)S(y, x)] + Tr[S(x, x)]Tr[S(y, y)] (1-16) 



For both u-quarks propagating in the zero mode we get, due to trace cyclicity: 




Tr[MxW(x)]Tr[My)*o(y)]} = o 



enforcing thus the Pauli principle. 
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Chapter 2 



T) c decay 

2.1 Introduction 

The charmonium system has played an important role in shaping our knowl- 
edge of perturbative and non-perturbative QCD. The discovery of the J/ip 
as a narrow resonance in e + e~ annihilation confirmed the existence of a new 
quantum number, charm. The analysis of charmonium decays in e + e~ pairs, 
photons and hadrons established the hypothesis that the J/tp and t] c are, to a 
good approximation, non-relativistic 3 Si and 1 5'o bound states of heavy charm 
and anti-charm quarks. However, non-perturbative dynamics does play an 
important role in the charmonium system [9,10]. For example, an analysis of 
the ip spectrum lead to the first determination of the gluon condensate. 

The total width of charmonium is dominated by short distance physics and 
can be studied in perturbative QCD [11]. The only non-perturbative input in 
these calculations is the wave function at the origin. A systematic frame- 
work for these calculations is provided by the non-relativistic QCD (NRQCD) 
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factorization method [12]. NRQCD facilitates higher order calculations and re- 
lates the decays of states with different quantum numbers. QCD factorization 
can also be applied to transitions of the type ip' — > ip + X [13, 14]. 

The study of exclusive decays of charmonium into light hadrons is much 
more complicated and very little work in this direction has been done. Per- 
turbative QCD implies some helicity selection rules, for example r] c -j* pp,pp 
and J/ip -/> pn,pa\ [15,16], but these rules are strongly violated [17]. The 
J/ip decays mostly into an odd number of Goldstone bosons. The average 
multiplicity is ~ (5 — 7), which is consistent with the average multiplicity in 
e + e~ annihilation away from the J/ip peak. Many decay channels have been 
observed, but none of them stand out. Consequently, we would expect the 
rj c to decay mostly into an even number of pions with similar multiplicity. 
However, the measured decay rates are not in accordance with this expecta- 
tion. The three main decay channels of the r\ c are KKn, rjTnr and tj'tttt, each 
with an unusually large branching ratio of ~ 5%. Bjorken observed that these 
three decays correspond to a quark vertex of the form (cc)(ss)(dd)(uu) and 
suggested that rj c decays are a "smoking gun" for instanton effects in heavy 
quark decays [18]. 

We shall try to follow up on this idea by performing a more quantitative 
estimate of the instanton contribution to r\ c and Xc decays. In section 12.21 we 
review the instanton induced effective lagrangian. In the following sections we 
apply the effective lagrangian to the decays of the scalar glueball, eta charm, 
and chi charm. We should note that this investigation should be seen as 
part of a larger effort to identify "direct" instanton contributions in hadronic 
reactions, such as deep inelastic scattering, the A J =1/2 rule, or 77 production 



37 



in pp scattering [19-22]. 



2.2 Effective Lagrangians 

Instanton effects in hadronic physics have been studied extensively [23,24]. 
Instantons play an important role in understanding the U(1)a anomaly and 
the mass of the rj'. In addition to that, there is also evidence that instantons 
provide the mechanism for chiral symmetry breaking and play an important 
role in determining the structure of light hadrons. All of these phenomena are 
intimately related to the presence of chiral zero modes in the spectrum of the 
Dirac operator in the background field of an instanton. The situation in heavy 
quark systems is quite different. Fermionic zero modes are not important and 
the instanton contribution to the heavy quark potential is small [25]. 

This does not imply that instanton effects are not relevant. The non- 
perturbative gluon condensate plays an important role in the charmonium 
system [9,10], and instantons contribute to the gluon condensate. In general, 
the charmonium system provides a laboratory for studying non-perturbative 
glue in QCD. The decay of a charmonium state below the DD threshold 
involves an intermediate gluonic state. Since the charmonium system is small, 
r cc ~ (vnic)' 1 < Aq^ d , the gluonic system is also expected to be small. For 
this reason charmonium decays have long been used for glueball searches. 

Since charmonium decays produce a small gluonic system we expect that 
the cc system mainly couples to instantons of size p ~ r c5 ~ (vm c )~ l . In 
this limit the instanton effects can be summarized in terms of an effective 
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lagrangian 11.101 discussed in chapter [TJ 




(2.1) 



Expanding the effective lagrangian in powers of the external gluon field 
gives the leading instanton contribution to different physical matrix elements. 
If the instanton size is very small, p <C tti' 1 , we can treat the charm quark 



contributes to charmonium decays. Since the density of instantons grows as a 
large power of p the contribution from this regime is very small. In the realistic 
case p ~ (vrric)^ 1 we treat the charm quark as heavy and the charm contribu- 
tion to the fermion determinant is absorbed in the instanton density d(p) . The 
dominant contribution to charmonium decays then arises from expanding the 
gluonic part of the effective lagrangian to second order in the field strength 
tensor. This provides effective vertices of the form (GG)(u'y 5 u)(d'y 5 d)(s'y 5 s), 
(G 2 )(M7 5 M)((i75(i)(ss), etc. 

We observe that the Nf = 3 fermionic lagrangian combined with the glu- 
onic term expanded to second order in the field strength involves an integral 
over the color orientation of the instanton which is of the form J dUlUijU^) 5 . 
This integral gives (5!) 2 terms. A more manageable result is obtained by using 
the vacuum dominance approximation. We assume that the coupling of the 
initial charmonium or glueball state to the instanton proceeds via a matrix 



mass as light and there is an effective vertex of the form (uu)(dd)(ss)(cc) which 
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element of the form (0 ++ |G 2 |0) or (0 + \GG\0). In this case we can use 

<0 ++ |G^G>> = 12{N \_ 1 f b (V uP - 5^5 ua )(0 ++ \G a ; a G^\0) (2.2) 

in order to simplify the color average. The vacuum dominance approximation 
implies that the color average of the fermionic and gluonic parts of the inter- 
action can be performed independently. In the limit of massless quarks the 
instanton (7) and anti-instanton (A) lagrangian responsible for the decay of 
scalar and pseudoscalar charmonium decays is given by 

Ci+a = J dz ^ d P ( N f/ l)a { {° 2 ~ GG ) x L v + ( g2 + GG ) x Am} • 

(2.3) 

Here, CfjA is the color averaged Nf = 3 fermionic effective lagrangian [29,23, 
24]. 

2.3 Scalar glueball decays 

Since the coupling of the charmonium state to the instanton proceeds via an 
intermediate gluonic system with the quantum numbers of scalar and pseu- 
doscalar glueballs it is natural to first consider direct instanton contributions 
to glueball decays. This problem is of course important in its own right. Ex- 
perimental glueball searches have to rely on identifying glueballs from their 
decay products. The successful identification of a glueball requires theoretical 
calculations of glueball mixing and decay properties. In the following we com- 
pute the direct instanton contribution to the decay of the scalar ++ glueball 
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state into nir, KK, rjrj and rjr)'. 

Since the initial state is parity even only the G 2 term in equ. (|2.H|) con- 
tributes. The relevant effective interaction is given by 



C 



I+A 



dz d Q (p) 



dp 



p 5 iV c 2 - 1 V a, 
[(uu)(dd)(ss) + (u'j 5 u)(d'y 5 d)(ss) 

(M7 5 w)((ic/)(s7 5 s) + (uu)(d>y 5 d)(s>y 5 s)} 



+ - 



(ut a u)(dt a d)(ss) + (ut a -f 5 u)(dt a -f b d)(ss 



+ (ut a ^u)(dt a d)(s^s) + (ut a u)(dt a 7 5 ^)(s7 5 s) 
3 

--[(ut a a fW u)(dt a a^d)(ss) + (ut a a^ 5 u)(dt a a^ 5 d)(ss) 
+ (ut a a llul 5 u)(dt a a llu d)(s 1 5 s) + (ut a a^u)(dt a o-^ b d)(s~i b s)] 

-^d abc [(ut a a fiu u)(dt b a fiu d)(st c s) + (ut a a^u)(dt b a^d)(st c s) 
+ (ut a a^ 5 u)(dt b a^d)(st c -f 5 s) + {ut a a^u)(dt b a^d){st c 1 b s)] 
+ (2 cyclic permutations u ^ d <-> s) 
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(wt a «)(«)(st C s) + (Mt a 7 5 «)(rft 6 7 5 d)(st C S 



+ (Mt a 7 5 M)(^)(s 7 5 t c s) + (Mt a M)(dtVrf)(st c 7 5 s) 
9 



32 J 



uf'a^yu) {dt b a vl d) (st c a lljL s) 
+ {ut a a^u) (dt b a v ^d) (st c a^s) 
+ [ut a a^u) (dt b a ul d) {st c a lfl ^ 5 s) 
+ {ut a a^u) (dt b a v ^d) {st c a^s) 



(2.4) 



Let us start with the process ++ — ► 7T7t. In practice we have Fierz rearranged 
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equ. (|2.4|) into structures that involve the strange quark condensate ss as well 
as operators with the quantum numbers of two pions. In order to compute the 
coupling of these operators to the pions in the final state we have used PCAC 
relations 



(0|J 7 5 n|7r + ) = V '"* J * = K„ (2.5) 
m u + m d 

(0\s^u\K^) = iV * m " f « = K K . (2.6) 



The values of the decay constants are f n = 93 MeV, = 113 MeV [30]. We 
also use Q u = (uu) = —(248 MeV) 3 and Qd = Q u as well as Q s = 0.66Q U [31]. 
The coupling of the rj' meson is not governed by chiral symmetry. A recent 
analysis of rj — rj' mixing and the chiral anomaly gives [32] 

(0|t2 7 5u|77) = (0| d l5 d\ V ) = -z(358 MeV) 2 = K«, (2.7) 

(0\u^u\v') = (0| Hd\i) = -i(320 MeV) 2 = ATJ, (2.8) 

(0|s 75 s|77) = z(435 MeV) 2 = K^, (2.9) 

(0\sj 5 s\r)') = — i(481 MeV) 2 = K%. (2.10) 

Finally, we need the coupling of the glueball state to the gluonic current. This 
quantity has been estimated using QCD spectral sum rules [33, 34] and the 
instanton model [35]. We use 

(0 ++ \g 2 G 2 \0) = A = 15 GeV 3 . (2.11) 

We can now compute the matrix element for ++ — > n + 7T~ . The interaction 
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vertex is 

£#3~ = fdz /^)jv57i(^) Q^yxiCa.G 3 )^)^^)- 

(2.12) 

The integral over the position of the instanton leads to a momentum conserving 
delta function, while the vacuum dominance approximation allows us to write 
the amplitude in terms of the coupling constants introduced above. We find 

<0 ++ (<z)|7r+(p>-GO> = (2 7 r) 4 6 4 (q-p + -p-)-^XoQ s Kl (2.13) 

where 

-/>)^(f)(^v) s 

The instanton density do(p) is known accurately only in the limit of small p. 
For large p higher loop corrections and non-perturbative effects are important. 
The only source of information in this regime is lattice QCD [36-39]. A very 
rough caricature of the lattice results is provided by the parameterization 

^ = ln 5(p- Pc ), (2.15) 
p b 2 

with no ~ lfm -4 and p c ~ 0.33 fm. This parameterization gives a value of 
A = (379 MeV) -9 . Another way to compute A is to regularize the integral 
over the instanton size by replacing d(p) with d(p) exp(— ap 2 ). The parameter 
a can be adjusted in order to reproduce the size distribution measured on 
the lattice. We notice, however, that whereas the instanton density scales as 
p h ~ 5 ~ p 4 , the decay amplitude scales as p b+8 ~ p 17 . This implies that the 
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results are very sensitive to the density of large instantons. We note that when 
we study the decay of a small-size bound state the integral over p should be 
regularized by the overlap with the bound state wave function. We will come 
back to this problem in section I2~4l below. 

We begin by studying ratios of decay rates. These ratios are not sensitive 
to the instanton size distribution. The decay rate ++ — > tt + it~ is given by 



r +H 



1 sfm 



0+H 



Ami 



16tt 



A 

167T 



n 2 



(2.16) 



The decay amplitude for the process ++ — > 7i n is equal to the ++ — > 
7t + tt~ amplitude as required by isospin symmetry. Taking into account the 
indistinguishability of the two ttq we get the total nn width 
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ml ++ - Ami 



32n 
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167T 



(2.17) 



In a similar fashion we get the decay widths for the KK, rp~], rjr)' and rfrj' 
channels 
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Glueball decay rates vs glueball mass, ccrrespbnding tb p ave =0.29 fm and m s =0 MeV 
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Figure 2.1: Scalar glueball decay rates plotted a function of the mass of 
the scalar glueball. The rates shown in this figure were computed from the 
instanton vertex in the chiral limit. The average instanton size was taken to 
be p = 0.29 fm. 



Here, KK refers to the sum of the K + K~ and KqKq final states. We note that 
in the chiral limit the instanton vertices responsible for nir and KK decays are 
identical up to quark interchange. As a consequence, the ratio of the decay 
rates T ++^ 7T7T /T ++_ ¥K ^ is given by the phase space factor multiplied by the 
ratio of the coupling constants 



(2.18) 

The main uncertainty in this estimate comes from the value of m s , which is 
not very accurately known. We have used m s = (140 ± 20) MeV. The ratio of 
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txtx to r\r\ decay rates is not affected by this uncertainty, 



r ++_ 



0.69, ' 0+H 



4m^ 



4m^ 



(2.19) 



In Fig l2.ll we show the decay rates as functions of the glueball mass. We 
have used Aqcd — 300 MeV and adjusted the parameter a to give the average 
instanton size p = 0.29 fm. We observe that for glueball masses mo++ > 1 
GeV the KK phase space suppression quickly disappears and the total decay 
rate is dominated by the KK final state. We also note that for m ++ > 1.5 
GeV the rjrj rate dominates over the tttt rate. 

In deriving the effective instanton vertex equ. (|2.12|) we have taken all 
quarks to be massless. While this is a good approximation for the up and 
down quarks, this it is not necessarily the case for the strange quark. The 
m s ^ contribution to the effective interaction for ++ decay is given by 



c 



p 5 



3^4 



1 TX 6 p 



7r p m s p(a s G ) x 



If- 3 

- I (uu)(dd) + (u-f 5 u)(d-f 5 d) + - [(ut a u)(dt a d) + (u-fH a u)(d-fH a d) 

3-3 
- -(ua^ v t a u)(da^t a d) - -(ua^h a u)(da^H a d) 



There is no m s ^ contribution to the KK channel. The m s ^ correction 
to the other decay channels is 
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Figure 2.2: Same as Fig. I2.ll but with m s ^ corrections in the instanton 
vertex taken into account. The results shown in this figure correspond to 
m s = 140 MeV. 
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where 



B 



%k{p> 



a 2 



-n 2 p 3 l p . 



(2.20) 



N 2 -l\a 2 s J\3 
The decay rates with the m s ^ correction to the instanton vertex taken into 
account are plotted in Fig. 12.21 We observe that effects due to the finite strange 
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Figure 2.3: Dependence of glueball decay rates on the average instanton size. 
The results shown in this figure correspond to the instanton vertex with m s ^ 
terms included. The strange quark mass was taken to be m s = 140 MeV. 

quark mass are not negligible. We find that the irir , 7777', and 7/7/ channels 
are enhanced whereas the r\r\ channel is reduced. For a typical glueball mass 
m o++ = (1-5 — 1.7) GeV the ratio r = B(tttt) / 'B(KK) changes from r ~ 0.25 
in the case m s = to r ~ 0.55 for m s ^ 0. In Fig. l2.3l we show the dependence 
of the decay rates on the average instanton size p. We observe that using the 
phenomenological value p = 0.3 fm gives a total width r ++ — 100 MeV. We 
note, however, that the decay rates are very sensitive to the value of p. As 
a consequence, we cannot reliably predict the total decay rate. On the other 
hand, the ratio of the decay widths for different final states does not depend on 
p and provides a sensitive test for the importance of direct instanton effects. 

In Tab. l2.ll we show the masses and decay widths of scalar-isoscalar mesons 
in the (1-2) GeV mass range. These states are presumably mixtures of mesons 
and glueballs. This means that our results cannot be directly compared to 
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resonance 


full width r(MeV) 


Mass (MeV) 


decay channels 


/o(1370) 


200-500 


1200-1500 


pp dominant 
7T7T, KK, rjrj seen 


/o(1500) 


109 ±7 


1507 ±5 


V T KK = 0.19 ±0.07 
^ = 0.095 ±0.026 
J? = 0.18 ±0.03 


/o(1710) 


125 ± 10 


1713 ±6 


= 0.39 ±0.14 
?m = 0.48 ±0.15 

1 KK 



Table 2.1: Masses, decay widths, and decay channels for scalar-isoscalar 
mesons with masses in the (1 — 2) GeV range. The data were taken from [30]. 



experiment without taking into account mixing effects. It will be interesting 
to study this problem in the context of the instanton model, but such a study 
is beyond the scope of this work. It is nevertheless intriguing that the /o(1710) 
decays mostly into KK. Indeed, a number of authors have suggested that the 
/o(1710) has a large glueball admixture [40-43]. 



2.4 Eta charm decays 

The r] c is a pseudoscalar J PC = h charmonium bound state with a mass 
m Vc = (2979 ± 1.8) MeV. The total decay width of the rj c is T Vc = (16 ± 3) 
MeV. In perturbation theory the total width is given by 

rfo - 2g) = 3^ n (l ± 4.4-i) . (2.21) 

Here, ^(0) is the 1 So ground state wave function at the origin. Using m c = 1.25 
GeV and a s (m c ) = 0.25 we get 1^(0)1 — 0.19GeV 3//2 , which is consistent with 
the expectation from phenomenological potential models. Exclusive decays 
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cannot be reliably computed in perturbative QCD. As discussed in the intro- 
duction Bjorken pointed out that r\ c decays into three pseudoscalar Goldstone 
bosons suggest that instanton effects are important [18]. The relevant de- 
cay channels and branching ratios are B(KKir) = (5.5 ± 1.7)%, B(t]7tit) = 
(4.9 ± 1.8)% and B^'nir) = (4.1 ± 1.7%). These three branching ratios are 
anomalously large for a single exclusive channel, especially given the small 
multiplicity. The total decay rate into these three channels is (14.5 ± 5.2)% 
which is still a small fraction of the total width. This implies that the as- 
sumption that the three-Goldstone bosons channels are instanton dominated 
is consistent with our expectation that the total width is given by pertur- 
bation theory. For comparison, the next most important decay channels are 
B{2(ix + ix~)) = (1.2 ±0.4)% and B(pp) = (2.6 ±0.9)%. These channels do not 
receive direct instanton contributions. 

The calculation proceeds along the same lines as the glueball decay calcula- 
tion. Since the t] c is a pseudoscalar only the GG term in equ. ()2.Hj) contributes. 
The relevant interaction is 




) 



3 




X 



+ (uu)(dd)(s'j 5 s) + (wy 5 u)(d'j 5 d)(s'y 5 s)} 

+- {ut a ^u){dt a d){ss) + (ut a u){dt a ^d){ss) + (ut a u)(dt a d)(s^s) 



+ {ut a ^u){dt a ^d){s^s) 
3 

--[(ut a a^ 5 u)(dt a a^d)(ss) + (ut a a^u)(dt a a^ 5 d)(ss) 



+ (ut a a^u){dt a a^d)(srfs) + (ut a a^u)(dt a a^d)(s^ s)] 
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X 



^d abc [{ut a a^u){dt b a^d){st c s) + (ut a a liu u)(dt b a liul 5 d)(st c s) 
+ (ut a a^u)(dt b a fMU d)(st c 1 5 s) + (ut a a^ 5 u)(dt b a^ 5 d)(st c 1 5 s)] 
+ (2 cyclic permutations u <-> d <-> s) 

(ut a ^u)(dt b d){st c s) + (Mt a M)(Jt b 7 5 ^)(^ c s) 

+ (ut a u)(^)(s 7 5 t c s) + (ut a ^u)(dt b ^d)(st c ^ s) 

( , ut a (T At[/ 7 5 M)((it f> (Tj y7 (i)(st c (T 7M s) + (ut a a^ L ,u)(dt b a„ 1 "f 5 d)(st c a lfl s) 



■—d abc 
40 



q 

i f a* x 

32 



+ (ut a a^ w u)(dt b a u ^d)(st c a lfJi 'j 5 s) 
+ {ut a a^ 5 u)(dt b a u ^ 5 d)(st c a lfl r 



(2.22) 



The strategy is the same as in the glueball case. We Fierz- rearrange the la- 
grangian ()2.22|) and apply the vacuum dominance and PCAC approximations. 
The coupling of the rj c bound state to the instanton involves the matrix element 



A* = (vc\g 2 GG\0). 



(2.23) 



We can get an estimate of this matrix element using a simple two-state mixing 
scheme for the r\ c and pseudoscalar glueball. We write 



|?7 C ) = cos(0)|cc) +sin(0)\gg), 
|(T + ) = -sin(0)|cc) + cos(0)|0s). 



(2.24) 
(2.25) 



The matrix element f ric = {Q\2m c cy^c\ri c ) ~ 2.8 GeV 3 is related to the char- 
monium wave function at the origin. The coupling of the topological charge 
density to the pseudoscalar glueball was estimated using QCD spectral sum 
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rules, A -+ = (0\g 2 GG\0~ + ) ~ 22.5 GeV 3 [34]. Using the two-state mixing 
scheme the two "off-diagonal" matrix elements /q-+ = (0|2m c c75c|0 h ) and 
A r?c = (0\g 2 GG\rj c ) are given in terms of one mixing angle 9. We can estimate 
this mixing angle by computing the charm content of the pseudoscalar glueball 
using the heavy quark expansion. Using [44] 

C7 5 c = ^—a s GG + O ( , (2.26) 
o7rm c \ m cJ 

we get /o-+ — 0.14 GeV 3 and a mixing angle 9 ~ 3°. This mixing angle 
corresponds to 

X Vc ~ 1.12 GeV 3 . (2.27) 

The uncertainty in this estimate is hard to assess. Below we will discuss a 
perturbative estimate of the instanton coupling to t] c . In order to check the 
phenomenological consistency of the estimate equ. ()2.27|) we have computed 
the 7] c contribution to the (g 2 GG(0)g 2 GG(x)) correlation function. The results 
are shown in Fig. 12.41 The contribution of the pseudoscalar glueball is deter- 
mined by the coupling constant A -+ introduced above. The couplings of the rj, 
t]' and 7/(1440) resonances can be extracted from the decays J/ip — > 777 [45] . We 
observe that the r\ c contribution is strongly suppressed, as one would expect. 
We also show the t] c and h glueball contributions to the (0750(0)0750(2;)) cor- 
relation function. We observe that even with the small mixing matrix elements 
obtained from equs. (j2.24H2.25j) the glueball contribution starts to dominate 
the r] c correlator for x > 1 fm. 

We now proceed to the calculation of the exclusive decay rates. There 
are four final states that contribute to the KKn channel, r\ c — > K + K~n°, 
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x [fm] 



Figure 2.4: Resonance contributions to the pseudoscalar glueball correlation 
function {g 2 GG(0)g 2 GG(x)) and the charmonium correlator (0750(0)0750(0;)). 
Both correlation functions are normalized to free field behavior. In the case of 
the gluonic correlation function we show the glueball contribution compared 
to the 77, 7/, ?7(1440) and r\ c contribution. For the charmonium correlation 
function we show the rj c and glueball contribution. 
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K°K°7T , K + K°ir~ and K~K°ir + . Using isospin symmetry it is sufficient to 
calculate only one of the amplitudes. Fierz rearranging equ. (|2.22jl we get the 
interaction responsible for the rj c — > K + K~tt° 



C 



K+K-jt 
I+A 



dz 



p 5 



iV c 2 -l 



a? 



'-7r 2 p 3 



x - (a s GG) (s 7 5 u) (u 7 5 s) (d^d) . 



(2.28) 



The decay rate is given by 
^k+k-tt° = / (phase space) x |M| 2 



A\ c K^K 2 K x(O.lllMeV), 

(2.29) 

with A given in equ. ([2. 14)1 . Isospin symmetry implies that the other KKn 
decay rates are given by 



K+K-TT — 1 K°K°w° 



"^/Tj ) (2.30) 



The total KKn decay rate is 
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6 x 



x (0.111 MeV) 



(2.31) 



In a similar fashion we obtain 
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x (0.135 MeV), (2.32) 

2 

x (0.0893 MeV), (2.33) 

2 

x (0.0788 MeV), (2.34) 
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r\ G main decay rates vs average inslanton size 
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Figure 2.5: Decay widths t] c — > KKtx and rj c — > rjirir as a function of the 
average instanton size p. The short dashed line shows the experimental KKir 
width. 



KKrf 

r 

1 r]T]r] 





1 


2 x 






16-7T 


1 


" 3! 


- X 




6 


16tt 



-A\ c K^K\ 
-AX Vc (KffK s v 
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(2.35) 



n 2 



x (0.0698 MeV). (2.36) 



Here, the first factor is the product of the isospin and final state symmetriza- 
tion factors. The second factor is the amplitude and the third factor is the 
phase-space integral. 

In Fig. 12.51 we show the dependence of the decay rates on the average 
instanton size. We observe that the experimental KKir rate is reproduced for 
p = 0.29 fm. This number is consistent with the phenomenological instanton 
size. However, given the strong dependence on the average instanton size it is 
clear that we cannot reliably predict the decay rate. On the other hand, the 
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following ratios are independent of the average instanton size 
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riTTTT 
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where we have only quoted the error due to the uncertainty in m s . These 
numbers should be compared to the experimental results 



KKtt 



rinir 



rjTVTV 



exp 



exp 



1.1 ±0.5 

1.2 ±0.6. 



(2.42) 
(2.43) 



We note that the ratio B(j]inr)/ B^'-irir) is compatible with our results while 
the ratio B(KKn) / Bijinn) is not. This implies that either there are con- 
tributions other than instantons, or that the PCAC estimate of the ratio of 
coupling constants is not reliable, or that the experimental result is not re- 
liable. The branching ratios for rjTTTT and ri'itir come from MARK II/III ex- 
periments [46,47]. We observe that our results for B{KKr\)jB(KK'K) and 
B(KKrj') I 'B(KKtt) are consistent with the experimental bounds. 

Another possibility is that there is a significant contribution from a scalar 
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resonance that decays into txtx. Indeed, instantons couple strongly to the 
cr(600) resonance, and this state is not resolved in the experiments. We have 
therefore studied the direct instanton contribution to the decay t] c — > at]. After 
Fierz rearrangement we get the effective vertex 



C ar) = I dA (a s GG) ^ [(u'j 5 u)(dd)(ss) + (uu)(d-f 5 d)(ss) + (uu)(dd)(s-f 5 s)] 
dB m s (a s GG) - [{u^u){dd) + {uu)(d^d)] , (2.44) 



where the integrals A and B are defined in equ. (j2.14l2.20J) . The only new 
matrix element we need is f a = (a\uu + d<i|0) ~ (500 MeV) 2 [48]. We get 
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X 



(2.45) 



fa\ c [{AQ s - 2Bm s )Kf 1 + AK°Q d ] 



16tt 



Compared to the direct decay r] c —* rjinr the r] c — > r/a channel is suppressed 
by a factor ~ (2tt 2 /mlj ■ (Q q f a /K 2 ) 2 ~ 1/100. Here, the first factor is due to 
the difference between two and three-body phase space and the second factor 
is the ratio of matrix elements. We conclude that the direct production of 
a a resonance from the instanton does not give a significant contribution to 
Vc — > r)(r)')7nr. This leaves the possibility that the ttw channel is enhanced by 
final state interactions. 

Finally, we present a perturbative estimate of the coupling of the t] c to 
the instanton. We follow the method used by Anselmino and Forte in order 
to estimate the instanton contribution to rj c — > pp [49]. The idea is that 
the charmonium state annihilates into two gluons which are absorbed by the 
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Figure 2.6: The Feynman diagram corresponding to the perturbative treat- 
ment of charmonium decay. 

instanton. The Feynman diagram for the process is shown in Fig J2.6l The 
amplitude is given by 




ufa) A^ikJA^ih), (2.46) 

where u(p) and v(p) are free particle charm quark spinors and A^ cl (k) is the 
Fourier transform of the instanton gauge potential 

A a ' d (k) = -i^f^m, m = 4(1- l -K 2 {kp){k P ) 

The amplitude for the charmonium state to couple to an instanton is obtained 
by folding equ. (j2.46J) with the r] c wave function ip{p). In the non-relativistic 
limit the amplitude only depends on the wave function at the origin. 

The perturbative estimate of the transition rate is easily incorporated into 
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the results obtained above by replacing the product AX Vc in equs. 
according to 



AX Vc - / ^d (p) ^-ttVJ (4tt) — |-|^(0)|/, c (p) x (2.4? 



with 



fc 2 $(A;)$(A; - 2 Pc 



I ( p ) = / d^k —— \^ £p zr. (2.49) 

Here, p c = (m c , 0) ~ (M^/2,0) is the momentum of the charm quark in the 
charmonium rest frame. We note that because of the non-perturbative nature 
of the instanton field higher order corrections to equ. ()2.48j) are only suppressed 
by g 2 (m^)/g 2 {p). 

The integral I Vc cannot be calculated analytically. We use the parameteri- 
zation 

IM K -ly'ff^y, ■ ■ (2.50) 
1 + B (m c p) 4 log(l + l/(m c p)) 

which incorporates the correct asymptotic behavior. We find that Aq = 0.213 
and i?o = 0.124 provides a good representation of the integral. In Fig. 12.71 we 
show the results for the r] c decay rates as a function of the average instanton 
size. We observe that the results are similar to the results obtained from the 
phenomenological estimate equ. (|2.27|) . The effective coupling (AX Vc ) differs 
from the estimate equ. ()2.27|) by about a factor of 3. The experimental KKn 
rate is reproduced for p = 0.31 fm. 
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il c perturbative and phenomenological decay rates vs average instanton size 
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Figure 2.7: Decay rates T(r) c — ► KKw) and T(rj c — > titttt) as a function of the 
average instanton size p. We show both the results from a phenomenological 
and a perturbative estimate of the cc coupling to the instanton. 

2.5 Chi charm decays 

Another interesting consistency check on our results is provided by the study 
of instanton induced decays of the \c into pairs of Goldstone bosons. The Xc 
is a scalar charmonium bound state with mass m Xc = 3415 MeV and width 
r Xc = 14.9 MeV. In a potential model the \ c corresponds to the 3 Pq state. 
In perturbation theory the total decay rate is dominated by cc — > 2g. The 
main exclusive decay channels are Xc — * 2(7r + 7r~) and Xc ~^ 7r + -n-~K + K~ with 
branching ratios (2.4 ± 0.6)% and (1.8 ± 0.6)%, respectively. It would be 
very interesting to know whether these final states are dominated by scalar 
resonances. We will concentrate on final states containing two pseudoscalar 
mesons. There are two channels with significant branching ratios, Xc — > tt + 7T~ 
and Xc -> K + K~ with branching ratios (5.0 ± 0.7) • 10 -3 and (5.9 ± 0.9) ■ 10~ 3 . 
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The calculation of these two decay rates proceeds along the same lines as 
the calculation of the ++ glueball decays. The only new ingredient is the Xc 
coupling to the gluon field strength G 2 . We observe that the total Xc decay 
rate implies that (0\2m c cc\x c ) = 3.1 GeV 3 ~ (0\2m c cij 5 c\r] c ) . This suggests 
that a rough estimate of the Xc coupling to G 2 is given by 



\ Xc = ( Xc \g 2 G 2 \0) ~ \ Vc = 1.12 GeV . 



(2.51) 



Using this result we can obtain the Xc decay rates by rescaling the scalar 
glueball decay rates equ. (j2.20H2.20j) according to 



Xc— >ml,m2 



r ++_ 



»ml,m2 



A 



Xc 



An 



(2.52) 



i 0++ ->m Xc 



where ml, m2 labels the two-meson final state. In Fig. l2.8l we show the depen- 
dence of the Xc decay rates on the average instanton size p. We observe that 
the experimental 7r + 7r~ decay rate is reproduced for p = 0.29 fm. In Fig. 12.91 
we plot the ratio of decay rates for 7r + 7r~ and K + K~ . Again, the experimental 
value is reproduced for p ~ 0.3 fm. 

Finally, we can also estimate the cc coupling to the instanton using the 
perturbative method introduced in section 12.41 In the case of the Xc we use 
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where R'(0) ~ 0.39GeV 5/2 is the derivative of the 3 P 

wave function at the 
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Chi decay rates in instanton-specific channels vs p. 





Figure 2.9: Ratio B(\c — * k+tt^)/ B(x c ~~ > K + K~) of decay rates as a function 
of the average instanton size. The dashed line shows the experimental value 
0.84. We also show the experimental uncertainty, as well as the uncertainty in 
the instanton prediction due to the the value of the strange quark mass. 
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X decay rates in instanton-specific channels vs p. 



for perturbative and phenomenelogical methpds 
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Figure 2.10: Decay rates T(x c 7r + 7r _ ) and r(x c — > K + K~) as a function 
of the average instanton size p. We show both the results from a phenomeno- 
logical and a perturbative estimate of the cc coupling to the instanton. 



origin and I Xc is the loop integral 



In Fig. 12.101 we compare the perturbative result with the phenomenological 
estimate. Again, the results are comparable. The experimental 7r + 7r~ rate is 
reproduced for p — 0.29 fm. 

2.6 Conclusions 

In summary we have studied the instanton contribution to the decay of a 
number of "gluon rich" states in the (1.5-3.5) GeV range, the scalar glueball, 
the t] c and the Xc- in the case of charmonium instanton induced decays are 
probably a small part of the total decay rate, but the final states are very 



$(&)$(|2p c - k\) 15(fc - Vcf + 3m;? + Ak 2 
k 4 (2 Pc - k) 4 (k-p c ) 2 + m 2 c 



(2.53) 
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distinctive. In the case of the scalar glueball classical fields play an important 
role in determining the structure of the bound state and instantons may well 
dominate the total decay rate. 

We have assumed that the gluonic system is small and that the instan- 
ton contribution to the decay can be described in terms of an effective local 
interaction. The meson coupling to the local operator was determined using 
PC AC. Using this method we find that the scalar glueball decay is dominated 
by the KK final state for glueball masses m ++ > 1 GeV. In the physically in- 
teresting mass range 1.5 GeV < m ++ < 1.75 GeV the branching ratios satisfy 
B(r)rj) : B(inr) : B(KK) = 1 : (3.3 ± 0.3) : (5.5 ± 0.5). 

Our main focus in this work are rj c decays into three pseudoscalar Gold- 
stone bosons. We find that the experimental decay rate T(?7 C — > KKtt) can be 
reproduced for an average instanton size p = 0.31, consistent with phenomeno- 
logical determinations and lattice results. This in itself is quite remarkable, 
since the phenomenological determination is based on properties of the QCD 
vacuum. 

The ratio of decay rates B^tj'tttt) : B^tttt) : B(KKtt) = 1:1: (4.2±1.3) is 
insensitive to the average instanton size. While the ratio B^tj'tttt) : B(riTcir) = 
1 : 1 is consistent with experiment, the ratio B(rj7rir) : BiKK) = 1 : (4.2±1.3) 
is at best marginally consistent with the experimental value 1.1 ±0.5. We have 
also studied Xc decays into two pseudoscalars. We find that the absolute decay 
rates can be reproduced for p = 0.29 fm. Instantons are compatible with the 
measured ratio B{K + K~) : ^tt+tt-) = 1.2 

There are many questions that remain to be answered. On the experimen- 
tal side it would be useful if additional data for the channels r\ c — > r]'im, r]im 
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were collected. One important question is whether (nir) resonances are im- 
portant. It should also be possible to identify the smaller decay channels 
r] c KKr], KKv( . In addition to that, it is interesting to study the distribu- 
tion of the final state mesons in all three-meson channels. Instantons predict 
that the production mechanism is completely isotropic and that the final state 
mesons are distributed according to three-body phase space. 

In addition to that, there are a number of important theoretical issues that 
remain to be resolved. In the limit in which the scalar glueball is light the decay 
++ — > TTir(KK) can be studied using effective lagrangians based on broken 
scale invariance [50-52]. Our calculation based on direct instanton effects is 
valid in the opposite limit. Nevertheless, the instanton liquid model respects 
Ward identities based on broken scale invariance [24] and one should be able 
to recover the low energy theorem. In the case ++ — > irn(KK) one should 
also be able to study the validity of the PCAC approximation in more detail. 
This could be done, for example, using numerical simulations of the instanton 
liquid. Finally we need to address the question how to properly compute 
the overlap of the initial cc system with the instanton. This, of course, is 
a more general problem that also affects calculations of electroweak baryon 
number violation in high energy pp collisions [53,54] and QCD multi-particle 
production in hadronic collisions [55]. 
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Chapter 3 

Instantons and the spin of the 
nucleon 

3.1 Introduction 

In this chapter we will try to understand whether the so-called "nucleon spin 
crises" can be related to instanton effects. The current interest in the spin 
structure of the nucleon dates from the 1987 discovery by the European Muon 
Collaboration that only about 30% of the spin of the proton is carried by the 
spin of the quarks [56] . This result is surprising from the point of view of the 
naive quark model, and it implies a large amount of OZI (Okubo-Zweig-Iizuka 
rule) violation in the flavor singlet axial vector channel. The axial vector 
couplings of the nucleon are related to polarized quark densities by 

g 3 A = Au — Ad, (3.1) 
g\ = Au + Ad-2As, (3.2) 
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= AS = Au + Ad + As. 



(3.3) 



The first linear combination is the well known axial vector coupling measured 
in neutron beta decay, g\ = 1.267 ±0.004. The hyperon decay constant is less 
well determined. A conservative estimate is g\ = 0.57 ± 0.06. Polarized deep 
inelastic scattering is sensitive to another linear combination of the polarized 
quark densities and provides a measurement of the flavor singlet axial coupling 
constant g A . Typical results are in the range g\ = (0.28 — 0.41), see [57] for a 
recent review. 

Since g A is related to the nucleon matrix element of the flavor singlet axial 
vector current many authors have speculated that the small value of g° A is in 
some way connected to the axial anomaly, see [58-60] for reviews. The axial 
anomaly relation 

PA} = ^G%G% + 2"W75<?/ (3.4) 

implies that matrix elements of the flavor singlet axial current A ^ are related 
to matrix elements of the topological charge density. The anomaly also im- 
plies that there is a mechanism for transferring polarization from quarks to 
gluons. In perturbation theory the nature of the anomalous contribution to 
the polarized quark distribution depends on the renormalization scheme. The 
first moment of the polarized quark density in the modified minimal subtrac- 
tion (MS) scheme is related to the first moment in the Adler-Bardeen (AB) 
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scheme by [61] 



AE^ = AE AB - Nf ^lAG(Q 2 ), (3.5) 

where AG is the polarized gluon density. Several authors have suggested 
that AT^ab is more naturally associated with the "constituent" quark spin 
contribution to the nucleon spin, and that the smallness of AS^j^ is due to a 
cancellation between AE^b an d AG. The disadvantage of this scheme is that 
AE^b is not associated with a gauge invariant local operator [62]. 

Non-perturbatively the anomaly implies that g° A = AS can be extracted 
from nucleon matrix elements of the pseudoscalar density mipi^^ip and the 
topological charge density g 2 G a ^ lu G a flu / '(327r 2 ) . The nucleon matrix element of 
the topological charge density is not known, but the matrix element of the 
scalar density g 2 G° ilv G a liV is fixed by the trace anomaly. We have [63] 

(N(p)\ J^G^W)) = C s (q 2 )m N u(p)u(p>), (3.6) 

with C s (0) = -1/b where 6 = 11- 2N f /3 is the first coefficient of the QCD 
beta function. Here, u(p) is a free nucleon spinor. Anselm suggested that 
in an instanton model of the QCD vacuum the gauge field is approximately 
self-dual, G 2 = ±GG, and the nucleon coupling constants of the scalar and 
pseudoscalar gluon density are expected to be equal, Gs(0) ~ Gp(0) [64], see 
also [65]. Using g° A = NfCp(Q) in the chiral limit we get g A ~ —Nf/b ~ —0.2, 
which is indeed quite small. 

A different suggestion was made by Narison, Shore, and Veneziano [66]. 
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Narison et al. argued that the smallness of AS = g\ is not related to the 
structure of the nucleon, but a consequence of the U(1)a anomaly and the 
structure of the QCD vacuum. Using certain assumptions about the nucleon- 
axial-vector current three-point function they derive a relation between the 
singlet and octet matrix elements, 



Here, f n = 93 MeV is the pion decay constant and x' to „(0) is the slope of the 
topological charge correlator 



with Qto P (x) = 9 2 G a llu G a pLl ,/ '(327T 2 ). In QCD with massless fermions topological 
charge is screened and Xto P (0) = 0. The slope of the topological charge corre- 
lator is proportional to the screening length. In QCD we expect the inverse 
screening length to be related to the if mass. Since the rf is heavy, the screen- 
ing length is short and x'topi®) * s small. Equation (|3.7|) relates the suppression 
of the flavor singlet axial charge to the large rf mass in QCD. 

Both of these suggestions are very interesting, but the status of the under- 
lying assumptions is somewhat unclear. In this work we would like to address 
the role of the anomaly in the nucleon spin problem, and the more general 
question of OZI violation in the flavor singlet axial-vector channel, by com- 
puting the axial charge of the nucleon and the axial- vector two-point function 
in the instanton model. There are several reasons why instantons are impor- 




(3.7) 




(3.8) 
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tant in the spin problem. First of all, instantons provide an explicit, weak 
coupling, realization of the anomaly relation equ. ()3.4j) and the phenomenon 
of topological charge screening [24,67]. Second, instantons provide a successful 
phenomenology of OZI violation in QCD [68]. Instantons explain, in particu- 
lar, why violations of the OZI rule in scalar meson channels are so much bigger 
than OZI violation in vector meson channels. And finally, the instanton liquid 
model gives a very successful description of baryon correlation functions and 
the mass of the nucleon [69, 70]. 

Our ideas are organized as follows. In Sect. 13.21 we review the calculation 
of the anomalous contribution to the axial-vector current in the field of an in- 
stanton. In Sect. 13 .31 and 13.41 we use this result in order to study OZI violation 
in the axial- vector correlation function and the axial coupling of a constituent 
quark. Our strategy is to compute the short distance behavior of the corre- 
lation functions in the single instanton approximation and to determine the 
large distance behavior using numerical simulations. In Sect. 13.51 we present 
numerical calculations of the axial couplings of the nucleon and in Sect. 13.61 we 
discuss our conclusions. Some results regarding the spectral representation of 
nucleon three-point functions are collected in appendix lAl 

3.2 Axial Charge Violation in the Field of an 
Instanton 

We would like to start by showing explicitly how the axial anomaly is realized in 
the field of an instanton. This discussion will be useful for the calculation of the 
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OZI violating part of the axial-vector correlation function and the axial charge 
of the nucleon. The flavor singlet axial-vector current in a gluon background 
is given by 

A ll (x) = Tr[ l5 ^S(x,x)} (3.9) 

where S(x, y) is the full quark propagator in the background field. The ex- 
pression on the right hand side of equ. (J3.9)) is singular and needs to be defined 
more carefully. We will employ a gauge invariant point-splitting regularization 



Tr [7 5 7 M 5'(x, x)] = hmTr 7 5 7 M S'(x + e,x- e)Pexp ^-i J A^dx) 

(3.10) 

In the following we will consider an (anti) instanton in singular gauge. The 
gauge potential of an instanton of size p and position z = is given by 



Here, f\ a is the 't Hooft symbol and R ab characterizes the color orientation 
of the instanton. The fermion propagator in a general gauge potential can be 
written as 

s(l , 9) = E !iMfiM, ,3,2) 

' A — m 

A 

where ^\(x) is a normalized eigenvector of the Dirac operator with eigenvalue 
A, Jp^\{x) = \^\(x). We will consider the limit of small quark masses. 
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Expanding equ. ()3.12|) in powers of m gives 

A^O A^O 

(3.13) 

Here we have explicitly isolated the zero mode propagator. The zero mode 
was found by 't Hooft and is given by 

*o(*) = - 7 2 ^ 1 2 ^3/ 2 7 ^ 7±0- (3.14) 

Here, (/> aa = e^/v^ is a constant spinor and 7± = (l±7 5 )/2 for an instanton/ 
anti-instanton. The second term in equ. ()3.13|) is the non-zero mode part of 
the propagator in the limit m — > [71] 

S» z (x,y) ee * X(X) * +Xiy) = p x A±(^vh± + A±(x,y) p y 7t (3.15) 

A^O 

where = — iA± and A±(x,y) is the propagator of a scalar field in the 
fundamental representation. Equ. ()3.15j) can be verified by checking that S± z 
satisfies the equation of motion and is orthogonal to the zero mode. The scalar 
propagator can be found explicitly 

A ±[x ,y) = *.(,,,,)-* (l + • V ) (3.16) 

where A = 1/ (47r 2 A 2 ) with A = x — y is the free scalar propagator. The 
explicit form of the non-zero mode propagator can be obtained by substituting 
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equ. (J3~TT)|) into equ. (j3~T5j) . We find 



1 + ^/1 + ^ 1 V ^ 

2 



A (z-!/) / P 



+ p2 + y 2 g: F ' X(T ± " A(T T " T^± " ^) } (3-17) 

Here, So = — Ay(27r 2 A 4 ) denotes the free quark propagator. As expected, the 
full non-zero mode propagator reduces to the free propagator at short distance. 
The linear mass term in equ. (j3.13|) can be written in terms of the non-zero 
mode propagator 



E * a( *].** (y) = / d*zS» z (x,z)S» z (z,y) = -A ± (x,y) 1± -A¥(x,yh T , 

A^O ^ 

(3.18) 

where A±(a;, y) is the scalar propagator and A±(x, y) = (x\(D 2 + a ■ G / '2)~ l \y) 
is the propagator of a scalar particle with a chromomagnetic moment. We 
will not need the explicit form of A±(x,y) in what follows. We are now in 
the position to compute the regularized axial current given in equ. ([3.10)1 . We 
observe that neither the free propagator nor the zero mode part will contribute. 
Expanding the non-zero mode propagator and the path ordered exponential 
in powers of e we find 



Tr [ l5 rS(x, x)] = ± 2( 2 (l\\ 3 > ( 3 - 19 ) 
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which shows that instantons act as sources and sinks for the flavor singlet axial 
current. We can now compare this result to the anomaly relation equ. (|H.4jl . 
The divergence of equ. (j3.19J) is given by 



The topological charge density in the field of an (anti) instanton is 

9 na no. _ _i_ jjfo (o o-]) 

16tt 2 7r 2 (x 2 + p 2 ) 4 ' 1 ; 

We observe that the divergence of the axial current given in equ. ()3.19|) does 
not agree with the topological charge density. The reason is that in the field 
of an instanton the second term in the anomaly relation, which is proportional 
to 771^75-0, receives a zero mode contribution and is enhanced by a factor 1/m. 
In the field of an (anti) instanton we find 



2m^ = T n2 ^f +p2Xi . (3.22) 



Taking into account both equ. (J3.21j) and (J3.22j) we find that the anomaly 
relation ()3.4j) is indeed satisfied. 
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3.3 OZI Violation in Axial- Vector Two-Point 
Functions 



In this section we wish to study OZI violation in the axial-vector channel due 
to instantons. We consider the correlation functions 

%,(x,y) = (j ll {x)j v (y)), (3.23) 

where is one of the currents 



v; = <H/r<v (p), v; = ^ ( w ), 

^ = h^r a ip (at), Al = ^7 M 75^ (/i), 



(3.24) 



where in the brackets we have indicated the mesons with the corresponding 
quantum numbers. We will work in the chiral limit m u = — > 0. The 
iso-vector correlation functions only receive contributions from connected di- 
agrams. The iso-vector vector (p) correlation function is 

(U s v )^(x,y) = 2(P£"V(x,j/) = -2(Tr [ ltl S(x, y)j u S(y, x)}). (3.25) 

The iso-singlet correlator receives additional, disconnected, contributions, see 
Fig. 13.11 The iso-singlet vector (u) correlator is given by 

(n° y V(x, y) = 2(prV(^ v) + <Pv s M*, V) (3.26) 
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Figure 3.1: Quark line diagrams that contribute to the vector and axial- vector 
two-point function in the iso-vector and iso-singlet channel. The solid lines 
denote quark propagators in a gluonic background field. The two diagrams 
show the connected and disconnected contribution. 

with 

(ifV(x,J/) = (Tr[^S(x,x)]Tr[ lv S(y,y)}). (3.27) 

The axial-vector correlation functions are defined analogously. At very short 
distance the correlation functions are dominated by the free quark contribution 
1T4 = 11^ = Uy = Uy ~ 1/x 6 . Perturbative corrections to the connected 
correlators are 0(ai s (x)/ir), but perturbative corrections to the disconnected 
correlators are very small, 0((a s (x)/ir) 2 ). In this section we will compute 
the instanton contribution to the correlation functions. At short distance, 
it is sufficient to consider a single instanton. For the connected correlation 
functions, this calculation was first performed by Andrei and Gross [72], see 
also [73]. Disconnected correlation function were first considered in [74] and a 
more recent study can be found in [75]. 

In order to make contact with our calculation of the vector and axial- vector 
three-point functions in the next section, we briefly review the calculation of 
Andrei and Gross, and then compute the disconnected contribution. Using 
the expansion in powers of the quark mass, equ. (j3.13j) . we can write 

(P v on )± (x, V) = (Pv n )o U + A± + £± (3-28) 
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with 



(P^X = -Tr[rS Q (x,y)YS (y,x)}, (3.29) 
Ag = -TrlrS^ z (x,y)YS^ z (y,x)]-(P^(x,y), (3.30) 
Sf = -2Tr[ 7 ^ 0± ( x )vl>+ ± ( y ) 7 -A ± ( 2/ ,x)7 ± ]. (3.31) 

Using the explicit expression for the propagators given in the previous section 
we find 



AT 



p 2 h x h y r g^^p 



± 2vr 4 A 4 I 



p h x h y (2A a A f3 - A 5 Q/3 ) + h y (ypA a + yaA^) 
/^(iA + x^)] =f 2e^(/i 2/ A a y /3 - ^A^)} (3.32) 



and 



= -^hlh 2 y {(x-y + p 2 )5^ - (yV - x»y») T e^x a y p ] (3.33) 

with h x = l/(x 2 + p 2 ), A = x - y, and = - g^g"? + ^w 3 ^. 

Our result agrees with [72] up to a color factor of 3/2, first noticed in [76], 
a '-' sign in front of the epsilon terms, which cancels after adding instantons 
and anti-instantons, and a '-' sign in front of the 2nd term in B^ v . This sign 
is important in order to have a conserved current, but it does not affect the 
trace pw. Summing up the contributions from instantons and anti-instantons 



2 
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we obtain 



ft 



-{h x h y ) 2 p 2 



2YS 



(2tt 2 ) 2 A 8 2vr 4 
2Y, ■ A 

H — — (E^A^ + A^E" - E ■ A^; 



A 4 



A^A 1 



A 4 



+— (A 2 g^ - A^ l A u - A^ u + A V W) 
A 2 



(3.34) 



This result has to be averaged over the position of the instanton. We find 



7l z 



2 



dA^dAy 



G{A\p) + 2G\A\p)g 



(3.35) 



±(x - z,y - z) = — 



d 2 



OA 2 



G{A\p) + 2G\A\p) 



9' 
(3.36) 



with 



G\A\p) 



dG{A\p) 



P 



2 P 2 n 



A 2 



(3.37) 



<9A 2 A 4 

and £ 2 = A 2 /(A 2 + 4p 2 ). The final result for the single instanton contribution 
to the connected part of the vector current correlation function is 



24 p 4 d 



tt 2 A 2 OA 2 \A 2 



lo s^ — 7 



(3.38) 
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Figure 3.2: Instanton contributions to the disconnected axial-vector corre- 
lation function. The left panel shows the single-instanton (non-zero mode) 
contribution. The right panel shows the instanton-anti-instanton (fermion 
zero mode) contribution. 



where we defined 



The computation of the connected part of the axial-vector correlator is very 
similar. Using equs. ()3. 3013.31)) we observe that the only difference is the sign 
in front of B^ v . We find 



5 ( P conw = (pconw _ 2(P% n )^ = -— \20A 2 G" + 56G'] , (3.40) 

TT 2 



with G' given in equ. ()3.37|) 

We now come to the disconnected part, see Fig. 13.21 In the vector channel 
the single instanton contribution to the disconnected correlator vanishes [74]. 
In the axial-vector channel we can use the result for Tr[j5'y fl S(x, x)] derived 
in the previous section. The correlation function is 



'pdisyu, y) _ Ap\x-zY{y-zY 
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Summing over instantons and anti-instantons and integrating over the center 
of the instanton gives 



p 1 



d 2 



2ir 2 dA.dA u 



(3.42) 



and 



i pdis^fifi 



4p 4 

IT 2 



A 



dA 2 ) J 



F(A,0- 



(3.43) 



We can now summarize the results in the vector singlet (u) and triplet (p), as 
well as axial- vector singlet (fi) and triplet (ai) channel. The result in the p 
and OJ channel is 



(iily» = (ir^ 



12 f 24 o A 

■^ + 2 dpn( P )- 2 ^F(A,0. 



(3.44) 



In the ai, fi channel we have 



(nir 



:n o r 



12 



7T 4 A ( 

12 

7rW 



+ / dpn(p) 
+ / dpn(p) 



7T 

2 



(20A 2 G" + 56G") 
(20A 2 G" + 56G") 



i^!(2F' + A 2 F") 



(3.45) 



(3.46) 



In order to obtain a numerical estimate of the instanton contribution we use a 
very simple model for the instanton size distribution, n(p) = Uq5{p — p), with 
p = 0.3 fm and Uq = 0.5 fm~ 4 . The results are shown in Fig. 13.31 

We observe that the OZI rule violating difference between the singlet and 
triplet axial-vector correlation functions is very small and repulsive. We can 
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Figure 3.3: correlation functions in the p, u, a\ and fx channel. All correlation 
functions are normalized to free field behavior. The data points show results 
from a numerical simulation of the random instanton liquid. The dashed lines 
show the single instanton approximation. 

also see this by studying the short distance behavior of the correlation func- 
tions. The non-singlet correlators satisfy 

(n 3 y r = (IW { 1 + j dpn(p) + ...}, (3.47) 
(nir = {VtiT {l - *>x* f dpn(p) + ...}. (3.48) 

As explained by Dubovikov and Smilga, this result can be understood in terms 
of the contribution of the dimension d = 4 operators (g 2 G 2 ) and (mqq) in the 
operator product expansion (OPE). The OZI violating contribution 

{u OZI r = (n r _ (n 3 )W = _ (n 3 ){; , j dp<p) (3 49) 

is of 0(x 6 ) and not singular at short distance. Our results show that it remains 
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small and repulsive even if x > p. We also note that the sign of the OZI- 
violating term at short distance is model independent. The quark propagator 
in euclidean space satisfies the Weingarten relation 



S(x,yy = -y 5 S(y,x)~? 5 . 



(3.50) 



This relation implies that Tr[S(x, 3^)7^ 75] is purely real. As a consequence we 
have 



Since the path integral measure in euclidean space is positive this inequality 
translates into an inequality for the correlation functions. In our convention 
the trace of the free correlation function is negative, and equ. ()3.5H) implies 
that the interaction is repulsive at short distance. The result is in agreement 
with the single instanton calculation. 

We can also study higher order corrections to the single instanton result. 
The two-instanton (anti-instanton) contributions is of the same form as the 
one- instanton result. An interesting contribution arises from instanton- anti- 
instanton pairs, see Fig. 13.21 This effect was studied in [77]. It was shown that 
the instanton-anti-instanton contribution to the disconnected meson channels 
can be described in terms of an effective lagrangian 



lim I Tr[S(x, x)j^ 5 ]Tr[S(y, y)~?^ 5 ] } > 0. 



(3.51) 



2G 



(V'7^75^) 2 



(3.52) 



with 




(3.53) 
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Figure 3.4: Correlation functions in the p, u, a\ and j\ channel. All correlation 
functions are normalized to free field behavior. The data points show results 
from unquenched simulations of the instanton liquid model. 

where n(pi, p?) is the tunneling rate for an instanton-anti-instanton pair and 
Tja is the matrix element of the Dirac between the two (approximate) zero 
modes. We note that this interaction is also repulsive, and that there is no 
contribution to the uj channel. 

Numerical results for the vector meson correlation functions are shown in 
Figs. 13.31 and 13.41 The correlation functions are obtained from Monte Carlo 
simulations of the instanton liquid as described in [78,79]. We observe that 
OZI violation in the vector channel is extremely small, both in quenched and 
unquenched simulations. The OZI violating contribution to the f\ channel is 
repulsive. In quenched simulations this contribution becomes sizable at large 
distance. Most likely this is due to mixing with an rj' ghost pole. We observe 
that the effect disappears in unquenched simulations. The pion contribution 
to the ai correlator is of course present in both quenched and unquenched 
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simulations. 

Experimentally we know that the p and u, as well as the a\ and f± meson, 
are indeed almost degenerate. Both iso-singlet states are slightly heavier than 
their iso-vector partners. To the best of our knowledge there has been only 
one attempt to measure OZI violating correlation functions in the vector and 
axial- vector channel on the lattice, see [80]. Isgur and Thacker concluded that 
OZI violation in both channels was too small to be reliably measurable in their 
simulation. 

3.4 Axial Vector Coupling of a Quark 

In this section we wish to study the iso-vector and iso-singlet axial coupling of 
a single quark. Our purpose is twofold. One reason is that the calculation of 
the axial- vector three-point function involving a single quark is much simpler 
than that of the nucleon, and that it is closely connected to the axial-vector 
two-point function studied in the previous section. The second, more impor- 
tant, reason is the success of the constituent quark model in describing many 
properties of the nucleon. It is clear that constituent quarks have an intrinsic 
structure, and that the axial decay constant of a constituent quark need not 
be close to one. Indeed, Weinberg argued that the axial coupling of a quark is 
(g\)Q — 0.8 [81]. Using this value of (g^Q together with the naive SU(6) wave 
function of the nucleon gives the nucleon axial coupling g\ = 0.8 • 5/3 ~ 1.3, 
which is a significant improvement over the naive quark model result 5/3. It is 
interesting to study whether, in a similar fashion, the suppression of the flavor 
singlet axial charge takes place on the level of a constituent quark. 
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Figure 3.5: Physical interpretation of the gauge invariant axial-vector three- 
point function of a quark in terms of a weak light-quark transition in a heavy- 
light Qq meson. 

In order to address this question we study three-point functions involving 
both singlet and triplet vector and axial-vector currents. The vector three- 
point function is 



The axial-vector function (H'Xqq) is defined analogously. We should note that 
equ. (|3.54|) is not gauge invariant. We can define a gauge invariant correlation 
function by including a gauge string. The gauge string can be interpreted as 
the propagator of a heavy anti-quark, see Fig. 13.51 This implies that the gauge 
invariant quark axial-vector three-point function is related to light quark weak 
transitions in heavy-light mesons. 

The spectral representation of vector and axial- vector three-point functions 
is studied in some detail in the appendix [XJ The main result is that in the 
limit that t/ 4 ^> z 4 ^> x 4 the ratio 



(Il a VQQ )f(x,z,y) = (q a (x)V«(z)f(y)). 



(3.54) 



Tr[(n^ Q ) 375 73] ; g_A 

Tr [( n V-QQ)474] 9V 



(3.55) 
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Figure 3.6: Quark line diagrams that contribute to the vector and axial- 
vector three-point function of a constituent quark. The solid lines denote 
quark propagators in a gluonic background field. The two diagrams show the 
connected and disconnected contribution. 

tends to the ratio of axial-vector and vector coupling constants, Qa/qv- We 
therefore define the following Dirac traces 

(U a VQQ r(x,z,y) = Tr[(U a VQQ rY], (3.56) 

(n a AQQ r(x,z, y ) = Tv[{n\ QQ y lb f\. (3.57) 

As in the case of the two-point function the iso-triplet correlator only re- 
ceives quark-line connected contributions, whereas the iso-singlet correlation 
function has a disconnected contribution, see Fig. 13.61 We find 

(W VQQ Y»{x,z,y) = (P^ Q r(x,z,y) = (Tr[S{x,z)>fS(z,y)>f\) (3.58) 

and 

(U° VQQ r(x,z,y) = (P^ Q r(x,z,y) - 2(P^ Q r(x,z,y) (3.59) 

with 

(Pfy Q T u (x,z,y) = (Tr[S(x,y)^]Tr[S(z,z)Y\), (3.60) 
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as well as the analogous result for the axial- vector correlator. 

In the following we compute the single instanton contribution to these 
correlation functions. We begin with the connected part. We again write 
the propagator in the field of the instanton as Szm + S^z + S m where Szm 
is the zero-mode term, S^z is the non-zero mode term, and S m is the mass 
correction. In the three-point correlation function we get contribution of the 
type SnzSnz, SzMS m and S m SzM 

( -peon \fiv _ r>^v | / tjHU p/ii/ \ 

\ r A/VQQ) — NZNZ ' °A/V \ r ZMm r mZMl 

= Tr[S NZ (x,z)rS NZ (z,y)Y] 

+c A/v Ti[-^ (x)^(z)^(-A ± (z,y) 1± Y] 
+c A/v Ti[(-A ± (x,z) 1± )^(-^ (z)^(yWl (3.61) 

where the only difference between the vector and axial-vector case is the sign 
of 'ZMm' and 'mZM' terms. We have ca/v — il fo r vector (axial vector) 
current insertions. The detailed evaluation of the traces is quite tedious and 
we relegate the results to the appendix iBl 

Our main goal is the calculation of the disconnected correlation function, 
which is related to OZI violation. In the single instanton approximation only 
the axial- vector correlator receives a non-zero disconnected contribution 

(Pf° Q r = TT[S(x,y) l5 Y}Tr[S(z,z) l5l % (3.62) 

see Fig. 13.71 We observe that the second trace is the axial- vector current in the 
field of an instanton, see equ. (|3.19jl . As for the first trace, it is easy to see that 
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Figure 3.7: Instanton contributions to the disconnected axial- vector three- 
point correlation function of a quark. The left panel shows the single-instanton 
(non-zero mode) contribution. The right panel shows the instanton-anti- 
instanton (fermion zero mode) contribution. 

neither the zero mode part of the propagator nor the part of Snz proportional 
to the free propagator can contribute. A straight-forward computation gives 

Tr[5(x,y)757 ] = =F 



vr 2 (x - y) 2 x 2 y 2 + ~p^l + ^ \p 2 + x 2 p 2 + y 2 

(3.63) 

Combined with equ. (J3.19j) we obtain 



Ms 2 P VsV(s - y y 



i^AQQ 



vr 4 (x - y) 2 x 2 y 2 (z 2 + p 2 f ^/( 1 + £i)( 1 + £ !) 
p 2 + x 2 p 2 + y 2 



(3.64) 



which has to be multiplied by a factor 2 in order to take into account both 
instantons and anti-instantons. 



C 0.8 




Figure 3.8: Axial and vector three-point functions of a quark as a function of 
the separation between the two quark sources. The correlation functions are 
normalized to free field behavior. The data points show results from numer- 
ical simulations of the instanton liquid and the dashed lines show the single 
instanton approximation. 

Results for the vector three-point function (nygg) 44 (r, r/2, 0) and axial- 
vector three-point function (U°J l qq) 33 (t,t/2,0) are shown in Fig. 13.81 We ob- 
serve that the vector and axial-vector correlation functions are very close to 
one another. We also note that the disconnected contribution adds to the con- 
nected part of the axial-vector three-point function. This can be understood 
from the short distance behavior of the correlation function. The disconnected 
part of the gauge invariant three-point function satisfies 

hm \(Uf QQ ) 3:i (x,z,y)} = ]im\Tr[S(x,x) l3l5 ]Tr[S(z, 2)7375] j > 0. (3.65) 

This expression is exactly equal to the short distance term in the disconnected 
fi meson correlation function. The short distance behavior of the connected 
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Figure 3.9: Axial and vector three-point functions of a quark as a function of 
the separation between the two quark sources. The data points show results 
from an unquenched instanton simulation. 

three-point function, on the other hand, is opposite in sign to the two-point 
function. This is related to the fact that the two-point function involves one 
propagator in the forward direction and one in the backward direction, whereas 
the three-point function involves two forward propagating quarks. A similar 
connection between the interaction in the f\ meson channel and the flavor 
singlet coupling of a constituent quark was found in a Nambu-Jona-Lasinio 
model [82,83]. It was observed, in particular, that an attractive coupling in 
the fi channel is needed in order to suppress the flavor singlet (g^g. 

The same general arguments apply to the short distance contribution from 
instanton-anti-instanton pairs. At long distance, on the other hand, we expect 
that IA pairs reduce the flavor singlet axial current correlation function. The 
idea can be understood from Fig. 13 .7[ see [58,84,85]. In an IA transition a 
left-handed valence up quark emits a right handed down quark which acts to 
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Figure 3.10: Ratio of axial-vector to vector correlation functions of a con- 
stituent quark calculated in the instanton liquid model. The open points show 
point-to-point correlation functions while the solid point is the zero momen- 
tum (point-to-plane) limit. The figure shows the iso- vector and iso-singlet 
correlation functions. 

shield its axial charge. We have studied this problem numerically, see Figs. 13.81 
13.101 We find that in quenched simulations the flavor singlet axial three-point 
function is significantly enhanced. This effect is analogous to what we observed 
in the fx channel and disappears in unquenched simulations. We have also 
studied the axial three-point function at zero three- momentum q — 0. This 
correlation function is directly related to the coupling constant, see App. El 
We find that the iso- vector coupling is smaller than one, {g\)Q — 0.9, in 
agreement with Weinberg's idea. The flavor singlet coupling, on the other 
hand, is close to one. We observe no suppression of the singlet charge of a 
constituent quark. We have also checked that this result remains unchanged 
in unquenched simulations. 
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3.5 Axial Structure of the Nucleon 



In this section we shall study the axial charge of the nucleon in the instanton 
model. We consider the same correlation functions as in the previous section, 
but with the quark field replaced by a nucleon current. The vector three-point 
function is given by 

m NN )f {x,v) = (rf(0)V?{y)ff{x)). (3.66) 

Here, rj a is a current with the quantum numbers of the nucleon. Three-quark 
currents with the nucleon quantum numbers were introduced by Ioffe [86] . He 
showed that there are two independent currents with no derivatives and the 
minimum number of quark fields that have positive parity and spin 1/2. In 
the case of the proton, the two currents are 

Vi = e a bc(u a C-f fl u b )-f 5 'y^d c , r} 2 = e abc (u a C(T llu u b )^a pLl/ d c . (3.67) 

It is sometimes useful to rewrite these currents in terms of scalar and pseudo- 
scalar diquark currents. We find 

Vl = 2{e abc (u a Cd b ) l5 u c -e abc (u a C l5 d b )u c }, (3.68) 
r)2 = 4{e abc (u a Cd b ) l5 u c + e abc (u a C l5 d b )u c }. (3.69) 

Instantons induce a strongly attractive interaction in the scalar diquark chan- 
nel e afec (ti fe C75(i c ) [69,87]. As a consequence, the nucleon mainly couples to the 
scalar diquark component of the Ioffe currents 771,2- This phenomenon was also 
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Figure 3.11: Quark line diagrams that contribute to the axial- vector three- 
point function of the proton. The solid lines denote quark propagators in 
a gluonic background field. The lines are connected in the same way that 
the Dirac indices of the propagators are contracted. The iso-vector and iso- 
singlet correlation functions correspond to g\ = —g A = 1 and g\ = g d A = 
1, respectively. The disconnected diagram only contributes to the iso-scalar 
three-point function. 

observed on the lattice [88]. This result is suggestive of a model of the spin 
structure that is quite different from the naive quark model. In this picture 
the nucleon consists of a tightly bound scalar-isoscalar diquark, loosely cou- 
pled to the third quark [89]. The quark-diquark model suggests that the spin 
and isospin of the nucleon are mostly carried by a single constituent quark, 
and that g A — g A - 

Nucleon correlation functions are defined by Ii^(x) = (^ a (0)^(a;)), where 
a, (3 are Dirac indices. The correlation function of the first Ioffe current is 

IT a/3 (x) = 2e abc e a , b , c , ((^ 5 Sf(0,x) lul5 ) Tr \^Sf(0, x)-y„C{S* (0, x)) T c\ ) 

(3.70) 

The vector and axial-vector three-point functions can be constructed in terms 
of vector and axial-vector insertions into the quark propagator, 

(T v X(x,y) = Sf(0,y) llx S b /(y,x), (3.71) 
{T*)f\x,y) = Sf(0,y) W S b /(y,x). (3.72) 
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The three-point function is given by all possible substitutions of equ. (|3.71|) 
and (J3.72)) into the two-point function. We have 



x(| g^hM^)f'^y)i^) aP Tr [ lp sf(o,x) lcr c(Su (o,*)fc~ 

+ 2fi (7 P 75^ C '(0, x) lal5 ) ap Tr [y p (r^(x, y) la C(S b a b ' (0, x)) T C 
(7p7 5 5f(0, x) 7 .7 5 ) Q/3 Tr [ lp Sf{^ *) 7 .C(Sf (0, x)) T C 

xTr [g^S d u d (y, y) + g^Sf (y, y)] J), (3.73) 

where the first term is the vector insertion into the d quark propagator in the 
proton, the second term is the insertion into the uu diquark, and the third 
term is the disconnected contribution, see Fig. 13.111 The vector charges of the 
quarks are denoted by g v . In the case of the iso- vector three-point function 
we have g v — 1, g v — — 1 and in the iso-scalar case g v — g v — 1. 

Vector and axial-vector three-point functions of the nucleon are shown in 
Figs. 13.121 and 13.131 In order to verify that the correlation functions are dom- 
inated by the nucleon pole contribution we have compared our results to the 
spectral representation discussed in the appendix |XJ see Fig. 13.121 The nu- 
cleon coupling constant was determined from the nucleon two-point function. 
The figure shows that we can describe the three-point functions using the phe- 
nomenological values of the vector and axial-vector coupling constants. We 
have also checked that the ratio of axial-vector and vector current three-point 
functions is independent of the nucleon interpolating field for x > 1 fm. The 
only exception is a pure pseudo-scalar diquark current, which has essentially 
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Figure 3.12: Vector, axial- vector three-point functions of the nucleon and nu- 
cleon two-point function calculated in the instanton liquid model. All correla- 
tion functions are normalized to free field behavior. The results are compared 
to a simple pole fit of the type discussed in the appendix. 

no overlap with the nucleon wave function. 

The main result is that the iso-vector axial-vector correlation function is 
larger than the vector correlator. The corresponding ratio is shown in Fig. 13. 131 
together with the ratio of q = correlation functions. We find that the 
iso-vector axial coupling constant is g\ = 1.28, in good agreement with the 
experimental value. We also observe that the ratio of point-to-point correlation 
functions is larger than this value. As explained in the appendix, this shows 
that the axial radius of the nucleon is smaller than the vector radius. Taking 
into account only the connected part of the correlation function we find a 
singlet coupling g° A = 0.79. The disconnected part is very small, g A (dis) = 
— (0.02 ± 0.02). Assuming that As ~ Au(dis) = Ad(dis) this implies that the 
OZI violating difference g\ — g A is small. This result does not change in going 
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Figure 3.13: Ratio of axial- vector to vector correlation functions of the nucleon 
calculated in the instanton liquid model. The open points show point-to-point 
correlation functions while the solid point is the zero momentum (point-to- 
plane) limit. The figure shows the iso- vector, connected iso-singlet, and full 
iso-singlet axial- vector correlation functions. 

from the quenched approximation to full QCD. 

We have also studied the dependence of the results on the average instanton 
size, see Fig. 13.141 We observe that there is a slight decrease in the iso-singlet 
coupling and a small increase in the iso-vector coupling as the instanton size 
is decreased. What is surprising is that the disconnected term changes sign 
between p = 0.3 fm and p = 0.35 fm. The small value g A (dis) = — (0.02±0.02) 
obtained above is related to the fact that the phenomenological value of the 
instanton size is close to the value where g A (dis) changes sign. However, even 
for p as small as 0.2 fm the disconnected contribution to the axial coupling 
g° A (dis) = — (0.05±0.02) is smaller in magnitude than phenomenology requires. 
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Figure 3.14: Axial coupling constants of the nucleon as a function of the 
instanton size p with the instanton density fixed at (N/V) = lfm -4 . We 
show the iso- vector, connected iso-singlet, and full iso-singlet axial coupling 
constant. 

3.6 Conclusions 

The main issue raised by the EMC measurement of the flavor singlet axial 
coupling is not so much why g° A is much smaller than one - except for the 
naive quark model there is no particular reason to expect g A to be close to 
one - but why the OZI violating observable g° A — g A is large. Motivated by 
this question we have studied the contribution of instantons to OZI violation 
in the axial-vector channel. We considered the fi — a\ meson splitting, the 
flavor singlet and triplet axial coupling of a constituent quark, and the axial 
coupling constant of the nucleon. We found that instantons provide a short 
distance contribution which is repulsive in the f\ meson channel and adds to 
the gauge invariant flavor singlet three-point function of a constituent quark. 
We showed that the sign of this term is fixed by positivity arguments. 
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We computed the axial coupling constants of the constituent quark and the 
nucleon using numerical simulations of the instanton liquid. We find that the 
iso-vector axial coupling constant of a constituent quark is (g A )Q = 0.9 and 
that of a nucleon is g A = 1.28, in good agreement with experiment. The result 
is also in qualitative agreement with the constituent quark model relation 
g\ = 5/3-(g\)Q. The flavor singlet coupling of quark is close to one, while that 
of a nucleon is suppressed, g\ = 0.77. However, this value is still significantly 
larger than the experimental value g A = (0.28 — 0.41). In addition to that, we 
find very little OZI violation, As ~ Au(dis) ~ —0.01. We observed, however, 
that larger values of the disconnected contribution can be obtained if the 
average instanton size is smaller than the phenomenological value of p ~ 1/3 
fm. 

There are many questions that remain to be addressed. In order to under- 
stand what is missing in our calculation it would clearly be useful to perform 
a systematic study of OZI violation in the axial-vector channel on the lattice. 
The main question is whether the small value of g A is a property of the nu- 
cleon, or whether large OZI violation is also seen in other channels. A study of 
the connected contributions to the axial coupling constant in cooled as well as 
quenched quantum QCD configurations was performed in [90]. These authors 
find g A (con) = Au(con) + Ad(con) ~ 0.6 in both cooled and full configura- 
tions. The disconnected term was computed by Dong et al. [91]. They find 
Au(dis) + Ad(dis) ~ -0.24. 

In the context of the instanton model it is important to study whether 
the results for g° A obtained from the axial-vector current three-point function 
are consistent with calculations of g A based on the matrix element of the 
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topological charge density GG [92-95]. It would also be useful to further 
clarify the connection of the instanton liquid model to soliton models of the 
nucleon [96]. In soliton models the spin of the nucleon is mainly due to the 
collective rotation of the pion cloud, and a small value for g° A is natural [97,98]. 
The natural parameter that can be used in order to study whether this picture 
is applicable is the number of colors, N c . Unfortunately, a direct calculations 
of nucleon properties for N c > 3 would be quite involved. Finally it would be 
useful to study axial form factors of the nucleon. It would be interesting to see 
whether there is a significant difference between the iso-vector and iso-singlet 
axial radius of the nucleon. A similar study of the vector form factors was 
recently presented in [99]. 
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Chapter 4 



Group integration 



4.1 General remarks 

Group integration became an important tool in physics with the development 
of lattice gauge theories [101]. Since then it enjoyed more attention in other 
well, like instantons [102, 103] or random matrix theory, e.g. [104]. 
In this chapter we elaborate on an idea of Michael Creutz and present an 
algorithm for computing integrals over elements of a compact group. 

A lot of progress has been made in the direction of solving the generating 
functional 



which enables one to express any other integral of a general function f(u, u') 
as: 



Different methods have been used to compute (|4.1|) . Direct integration over 




(4.1) 




(4.2) 
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parametrized SU{2) and SU(3) was given in [105], and polar decomposition of 
U(N) and integration over angular variables in [106]. Another approach was to 
construct and solve a partial differential equation which Z(J, Jt) satisfies due 
to its invariance w.r.t. left and right actions of the group [107,108]. The author 
of [109] expanded Z(J, J') in powers of invariants and computed numerically 
the coefficients up to some order using again a differential equation. Character 
expansion (for recent work see [111]) proved to be a powerful method, used for 
U(N) in [110]. More recently [104] a expansion of integrals over monomials 
in U (N) and O(N) matrix elements was developed, with the leading and next- 
to-leading order terms computed using a standard integral over C^ 2 with a 
modified Gaussian measure. 

The method we propose lacks the elegance of character expansion of gener- 
ating functional. Its simplicity, however, makes it a handy tool for evaluating 
group integrals over any compact group and any representation. Moreover, it 
is easily implement able as a computer algorithm once the group invariants are 
known. The same idea may also be applied to evaluating integrals of tensorial 
structure over any manifold with a measure possessing some kind of symmetry. 

The flow of ideas has the following structure: in section 14.21 we introduce 
the main properties of Haar measure. We state the algorithm in section l431 and 
exemplify it on fundamental and adjoint representations of SU(N) in section 
14.41 and 14.51 respectively. 
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4.2 Properties of Haar measure 

The algorithm being proposed relies heavily on the properties of the Haar 
measure 1 : 

For any compact group there exists a unique left-right invariant measure, 
normalized to unity, such that: 



where u, g are elements of the group and du is the Haar measure. 

The left-right invariance of the measure is the analog of translational in- 
variance of the R n integral and for finite groups represents the invariability of 
the number of elements w.r.t. multiplication by a fixed element of the group. 
The examples of group measure [101]: 




(4.3) 



(4.4) 



(4.5) 



for 17(1) = {e 1 ' 



n < (j) < 7r}: 




(4.6) 



for SU (2) = {a + ia • r | + a 



1}: 




(4.7) 



for a more detailed discussion see e.g. [101] and references therein. 
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4.3 The algorithm for group integration 

We will present here the main idea, mentioned before in [108] but, to our 
knowledge, never carried out completely. 

For the purpose of integration we will not need an exact expression for the 
measure, like f)4.fij) . (j4.7j) . We will only use the general properties of the Haar 
measure, making the algorithm applicable to any compact group. 

We are interested in the most general integral of a monomial in group 
elements: 

I l 1 ...l m ,ki...k m = f duu u u klk ...u kmlm , (4.8) 

'left' 'right' indices 

where u represents for SU(N) or u + for SO(N) or u^ 1 for a general compact 
group. 

In order to exemplify the idea, we will consider first the simple example of 

= Jdu Uij u kl . (4.9) 
The right invariance (|4.4|) implies: 

4,i = 9kk> iff gj'j ■ (4.10) 

This shows that l\' k has to be an invariant of the group in 'right' indices j, k. 
The same argument with left invariance of the measure compels l\' k to be 
separately an invariant in 'left' indices i, I as well. There is only one second 
rank invariant tensor with upper and lower index for SU(N): 8\. Therefore, 
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the integral ()4.9|) has to be of the form: 



= ai S l 5 3 ■ 



The constant a x is found easily by contracting with <5| and using UyVp 1 = S\ 
and the unit normalization (J4.3|) : 



So the integral then is: 



4 fc = = a* , (4.ii) 



which is nothing else but the orthogonality relation for fundamental represen- 
tation. The non-existence of one-index invariant leads directly to 

= yduuy = . (4.12) 

The same steps as in the above examples can be followed for any integral of 
the type (|4.8j) . Let us summarize them in the form of the following algorithm: 

• Generate all the independent invariant tensors with the index structure 
of (|4.8|) . The result of integration will be take the form: 



const x (left indices invariant) x (right indices invariant). 



all invariants 



Use index pairs exchange symmetry(e.g. <-> 2 2 j2) to identify the 
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tensors with the same numerical coefficients. 

• Compute the independent coefficients by multiplying ()4.8|) with the cor- 
responding invariant tensors (e.g. contracting indices with S tensor). 

Throughout the this work we deal with SU(N). However, the same philosophy 
is applicable to any compact group. 



4.4 Integration over fundamental representa- 



The invariant tensors depend on the group (e.g. vary with N in SU(N)) as well 
as representation. For fundamental representation of SU(N) we have two basic 
invariants, out of which the rest is constructed: 8\ and e ll "' JJV . The invariance 
can be seen easily from the unitarity and unimodularity. 
The result of a particular integral will depend on N, e.g. 2 



Despite of having different invariants for different N, there are many cases the 
result has the same form, but coefficients depend on N. An example is the 



tion of SU(N) 




(4.13) 



integral with the same number of u and it's, schematically 



I n = Jdu(u w)",on 



which we will focus in this paragraph. 



2 the coefficient of ^ is determined by multiplication with the invariant tensor itself: 
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Let us then follow the proposed algorithm and compute I n . First we gener- 
ate all left-right invariant terms. Using symmetries of the measure we identify 
the terms with the same numeric coefficients. Finally we multiply both sides 
with the very same set of invariants and solve the obtained system of linear 
equations for the unknown coefficients. 

So what are the possible invariants we can construct? The building blocks 
are Kronecker's 5 and e tensors. 

First of all, there are (n\) 2 terms of the form 

( 6 x . . . x 5) x ( 8 x . . . x 5 ) , (4.14) 

n X left indices nxright indices 

i.e. a § for each up-down pair both in left and right indices. 

Then, for n > N we can use e^...^ or e ll " lN to form terms like 

( 5 X . . 5 y X e[ ower indices x e Upper mdlCeS )i e f t X (. . .) right . (4-15) 
n-N 

Notice that e tensors must come in pairs, both with upper and lower indices, 
left or right (or both). The reason is obvious: 5's require both an upper and 
a lower index; if an e with lower indices is used, then the only invariant left 
for upper extra indices is another e. One could thus have pairs of e in left and 
right indices, having different structure for different N. Fortunately, this is not 
the CclSG, clS products of e's can be expanded in products of 5's. 

The important observation is, that any pair of e's in upper and lower indices 
can be written as a completely antisymmetric sum of products of Kronecker 5 
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tensors: 

e ill2 ..,J^= £ (-l) P < 0l) ---C (W • (4-16) 

permutations P 

Therefore, independently on N, the final result will only contain products 
of <5's, it will be of the form of (j4.14j) . 

Having noticed this, we only need to generate all the possible (S) n x (S) n 
combinations, identify common coefficients and compute them, (the actual val- 
ues of coefficients will depend on N). 

We will present the calculation up to n = 3, for the higher order inte- 
grals one would use a computer program to generate terms and compute the 
coefficients. 

4.4.1 Case n=2 

We will treat the integral 

4tS = j du u nn u klh u l2n u k ^ (4.17) 

in some detail in order to unveil the ideas that are useful in the more involving 
case of higher n. 

The integral will contain (2!) x (2!) terms: 

) ♦'<( X ) • (418 > 

where we represent 8 by a line joining the indices. 3 Each line has to join an 

3 In all the subsequent diagrams the places of indices are the same as in the first diagram, 
so we omit them for the sake of simplicity 
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upper left (right) index with a lower left (right) index. No connections between 
left and right indices are allowed. The indices in the above graphs have always 
the same position as in the first parenthesis of (|4.18|) 

The (iiji) <-> (22J2) and (&1/1) ^ (^2^2) symmetry of (|4.17|) allows us 
to identify the equal coefficients among fe 1 ...&4. Under (iiji) <-> (^2^2) the 
diagrams (1 and 4) as well as (2 and 3) interchange. Therefore b\ = 6 4 and 
&2 = h- 

It is easy to identify the diagrams that have the same coefficients by the 
topology of diagrams. Take pairs of 1 left and 1 right index, e.g. lower 
(iiji), («2j2) or upper (hk\), (£2^2)- What matters is how many of these pairs 
of indices are interconnected: in b\ and 64 terms one only has (1 upper pair - 1 
lower pair) 5 connections. These two diagrams are pair-exchange equivalent. 
The 62 diagram is topologically different: there is a (2 upper pairs - 2 lower 
pairs) ^-interconnection: both pairs of upper indices (lk) are connected to both 
pairs of lower (ij) indices. This diagram is pair-exchange equivalent to 63 and 
falls in the same topological set of diagrams. This easy way to identify dia- 
grams with the same coefficients will prove very useful in calculating integrals 
with higher n. 

Computing the coefficients is straight-forward: contract with 'right' 5^ to 
get on one hand, from (|4.11j) : 

jy u n 32 



108 



and, on the other hand (|4.18l) : 



(N b x + b 2 )5 l ^8f 2 +(Nb 2 + &0<5§a£*g . 



Matching the two sides we get the coefficients. The result then is: 



N 2 - 1 



dll Ui 1 j 1 Ui 2 j2 r U J k2h 



N{N 2 -1) 

— - — \5 h S kl S l2 5 k2 + 5 l2 5 k2 5 h 5 kl ] \5 h 5 k2 5 h 5 kl + 5 h 5 kl S h S k2 } 

jY"2 _ r»l 31 12 B 1 n Jl 13 32l N(N 2 — 1) 1 Jl 12 3 2 '131 12 bJ 



(4.19) 



4.4.2 Case n=3 



The number of terms in 



jhhh,kxk 2 h z Iduu- ■ u klh v- ■ v k2h v- ■ v k3h 



(4.20) 



increases to (3!) x (3!) = 36. They will be grouped in 3 topologically different 
classes having the same coefficients in the front: 

• 6 terms of type (1-1,1-1,1-1), like e.g. 

/ 1] I2 It2 h ki\ / \ 




'1 Jl '2 12 >3 13 



linking 1 pair of upper to 1 pair of lower indices. 



18 terms of type (2-2, 1-1) having 2 pairs of lower and 2 pairs of up- 
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per indices interconnected in addition to one upper to one lower pair 5 
connection. Examples of diagrams from this category: 



ij kj i 2 k 2 i 3 k 3 





H Jl '2 J2 '3 J 3 



12 terms of type(3-3) having all 3 pairs of upper indices and 3 pairs of 
lower indices interconnected, e.g. diagrams: 



li ki Lo k r\ l^k 





H Jl 1 2 J 2 '3 J 3 



The coefficients can be found in a similar way to the case of n — 2. The 
result is: 



du ^i^j/l ^ J2 ^ 33 ^ 



k2h 



Cl 



C 2 



+ c 3 



It] 1 2 k- 2 1 3 ^ 



' 1 x 2 J 2 '3 J3 



J 1 x 2 J 2 '3 J 3 



kj 1 ^ k. 2 1 k 



. . . altogether 6 terms 
. . . altogether 18 terms 
. . . altogether 12 terms 



+ 



+ 



J 1 '2 J 2 "3 J 3 



(4.21) 



where different terms in . . . are generated by interchanges of (Ik) pairs or 
(ij) pairs, and have the same topology as the first term in parenthesis. 
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For SU (N > 3) the coefficients are: 

(N 2 - 2) 

c 2 : 

° 3 = N(N 2 - l)(iV 2 -4) ^ 4 ' 22 ^ 

An important requirement is, that the invariant tensors are independent. 
The 3! terms of type S^S 1 ^ 5j* and permutations are indeed independent for 
SU(N), N > 3. However, for N = 2 there are 'too many' indices and only 2 
possible values for each index. This leads to the following relationship: 



N(N 2 


- 1)(N 2 
1 


-4) 


(N 2 


- 1)(N 2 
2 


-4) 


N(N 2 


- 1)(N 2 


-4) 



D = f ■ ■ ■ f lll2li = 

'l '2 '3 

(4.23) 



In this case of having dependent tensors one can follow the same strategy up to 
one point: after contracting the indices in order to figure out the coefficients, 
one has to express everything in terms of a basis of independent tensors, and 
only then compare the coefficients. For SU(2) this leads to the same form of 
result 1)4.21)1 with the parameter-dependent coefficients: 

1 

c i = t+ a 
o 

1 

Co = —t 

24 

c 3 = t (4.24) 

The parameter t cancels if one goes to the basis of 5 independent SSS tensors 
so we can as well set it to zero. 
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For higher n — 4, 5 . . . one has to take this subtlety into account also for 
SU(3), SU(4) . . . The same idea will prove useful in the next paragraph, when 
we treat the integral over adjoint representation. 



4.5 The integral over adjoint representation of 
SU(N) 

In a number of cases one encounters integrals of type 

V 2 . . . a n , bl b 2 • • A = j du ****** ' ' ' ^ • ^ 



left right 



where 



R ab = br r ^ X b uX a u) 



is the adjoint representation of SU(N) satisfying R ab X b = u\ a u, with A a the 
generators of SU(N) normalized to Tr(X a X b ) = 25%. 

To solve this integral, one can straight-forwardly integrate over u's and 
then contract the result with A's and perform the trace. For example: 



jduR ab = ^X b u X% jduu l3 u kl 

= hiX^SjA = ^TrX b TrX a = . (4.26) 

For higher n the number of terms of jdu (uu) n grows as (n!) 2 and the eva- 
luation of the integral, although straight-forward, becomes tedious. For this 
reason, we present here a more manageable way of computation, based on the 
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above algorithm. 

First of all, let's note some properties of R ab matrices: 

• R ab is the adjoint representation of SU(N) group 4 : 

R±> = K' Ru (4-27) 

• The matrices form a subgroup of SO(N 2 — 1): 

R ab R cb = R ab( R +yc = S ac (^g) 

• Under the left-right action of the group, it rotates as follows: 

R ab u—*gjih R ab R ac ^ R^ c R a< ^ R^ C R° b (4 29) 

The main idea of the algorithm applies here as well: the left-right invariance 
of the measure in ()4.25|) implies that 

l ai ...a n M...K = R-r' 1 ■ ■ ■ R a h n< R b g lb 'l ■ ■ ■ R b fi x l a ,.. <iK .. K . (4.30) 

In other words, I ai ...a n ,bi...b n h &s to be of the form 

const x (invariant in left indices) x (invariant in right indices) 

all invariants 

(4-31) 

4 We will denote by R^ b the adjoint matrix corresponding to the group element u. When 
there is no ambiguity on what element the R corresponds to, we will skip writing u's, usually 
it will mean that all the R matrices correspond to the same group element 
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Up to now, everything is similar to what we had for integration over funda- 
mental representation. However, the construction of invariants will be more 
involved, due to two extra 'building blocks': d abc and f abc . The SU(3) invari- 
ants up to the 5th rank can be found e.g. in [112, 114]. One might also find 
useful some SU(N) relations in [113]. 

The symmetries used for determining the common coefficients of invariants 

are: 

• (dibi) ( a jbj) exchange symmetry 

• a <-> b symmetry stemming from inversion invariance of the measure, as 



Following the algorithm presented above, we computed the following inte- 
grals: 

4.5.1 Case n=l 

The non-existence of 1 index invariant leads to 



in (jg~HD 




(4.32) 



as we explicitly computed in (j4.2fi|) . 
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4.5.2 Case n=2 

There is only one 5 second rank invariant: 5, and one easily gets the orthogo- 
nality relation for the adjoint representation: 



du R aibl R a2b2 



1 



N 2 - 1 



(4.33) 



4.5.3 Case n=3 



There are 2 third rank invariants: d and f and so 



du R aibl R a2b2 R a3b3 



1 



+ 



N(N 2 - 1) 
N 



aia 2 a 3 ^616263 



yaia 2 a 3 j 



(N 2 -A)(N 2 - 1) 



d a ^ a H blb2bz . (4.34) 



For SU(2) we have d aia2a3 = and one just disregards the second term in the 
result above. 

4.5.4 Case n=4 

We will use the following diagrammatic representation: 



5 ai 

CL 2 



daia 2 ada 3 a^a 



1 "3 
f 

'2 a 4 



da\a 2 aj aj 



3 we need separately an invariant for left and one for right indices 
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The bracketed terms will mean invariants both in left ('a') and right ('b') 
indices. The places of indices is always the same as in the following examples: 



/*1 ■•< >': ■'• 



A A A fo 3 
^axa-ia^aza^a "6 4 "6 2 



a 2 a 4 b 2 b 4? 



or 



\ ' / \ "2 a 4 b 2~ b 4l 

and so on. For SU(2) this is still a simple integral, since the only independent 
invariant tensors are again products of 5(The structure constants tensors re- 
duce to product of <5's as well, due to (jJUBJ)). The straight-forward application 
of the algorithm leads for SU(2) to: 



du R aibl R a2b2 R a3h R aib4 



2 

15 

/ 



"J "3 b l b 3 • 



»7 a 4 h b 4l 



+ 



+ 



+ 



+ 



(4.35) 



1 

30 



+ 



+ 



/ \ 



+ 



Following [114], one can construct 9 independent invariants of 4th rank for 
adjoint SU(N > 3) and 8 for SU(3). 
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In terms of diagrams, the n = 4 result for SU (N > 3) is: 



du R aibl R a2b2 R a3h R a4h = K u 

r / \ 



• Wy yW f- 1 



+ #4i 
+ #12 
+ # 42 
+ # 42 
+ # 45 
+ # fi 



+ 



+ 



+ 



7 



+ 



/ \ 



\ ' 1 ' 7" 



) + (' 



) + ( 
W 

/ \ 

+ ( 

V \ 
>: : + 



\ 

+ 

\ 



+ 



o 1 

o ! ! 

" . \ 7 

+ ;< I + I 

V'' -■/ \- 



+ i 

\ / 

\ /v ~\ /~ 

+ A + 

/ \ A — / v— 

\^N-\^ (■■- 
I + (a -j + (- 

+ h i + 

A; '.;\ / \ / l '.A 

+ 3 i y t i + 3 ! ) - i y + 

yA Ay ^ ^y \ i Ay 
+ , 



; + | 

/ \ 



A A 



where the coefficients for the case of SU(N > 3) are: 



+ 



+ 



+ 



+ 



+ 



-/ A 



(4.36) 



#ii 

#12 
#41 
#42 



N 4 - 6N 2 - 24 
N 4 (N 2 - 9)(N 2 - 1) 

N 2 - 12 
~N 4 {N 2 -9)(N 2 - 1) 

iV 2 - 12 
N 3 (N 2 — 9)(iV 2 - 1) 

iV 2 -6 
~jV3(AT2 _ 9 )(iv 2 - 1) 
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K 44 



3N 4 - 29N 2 + 48 



2jV 2 (iV 2 -l)(iV 2 -4)(iV 2 -9) 
iV 4 - 15iV 2 + 24 

~2iV 2 (iV 2 - 1){N 2 - 4)(iV 2 - 9) 
1 



(4.37) 



2(iV 2 -4)(iV 2 - 1; 



The subtlety for SU(3) relies on an extra, SU(3)-specific relation between d 
tensor components [113]. Diagrammatically: 



(4.38) 



This reduces the basis of independent tensors to 8. One has to keep this in 
mind when making contractions to get the coefficients. The system of equa- 
tions will then be under-determined and the parametric 6 solution for SU(3) 



is: 



91 

6480 ~ 
13 

~2160 



9 45 



-K 



1 



12 



-K 



15 



108 
1 

180 
1 

80 



+ K 42 
+ K 45 



(4.39) 



For higher n = 5, 6 . . . one needs to take into consideration the analogous 



Since K42, cancel when one expresses all in terms of independent basis, we can as 
well choose them to be zero. 
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relations of the equation (ICTJ) for SU(4), SU(5). 



4.6 Conclusions 

We have presented a 'down to earth' algorithm for evaluating group integrals. 
While lacking the elegance of more group-theoretic approaches (like character 
expansion), the strength of our method is in its simplicity, allowing one to 
reach the result in an easy and short way. We illustrated the method on 
examples of integrals over SU(N) fundamental and adjoint representation, 
but the algorithm is not restricted to these particular groups. While some of 
the results for J du(uu) n we derived have already been known, the results for 
J du(R ab ) 4 , to our knowledge, have not been published. 
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Chapter 5 



Summary and outlook 



We have studied instanton contributions to processes with unusual spin-flavor 
structure. The first part of thesis focused on the decay of scalar and pseu- 
doscalar charmonium and glueball. We found an internally consistent picture 
that agrees well with instanton phenomenology and lattice results. The decay 
rates of studied processes are reproduced by instanton computation with the 
average size of instanton p = 0.3 fm. The ratios of decay rates do not de- 
pend on the average size of instantons and agree, with one exception, with the 
experimental values. 

A better understanding of rj c decay would come with more experimental 
data. For example, it would be useful to compare our predictions with the 
data on smaller decay channels rj c — > KKrj, KKrj'. The distribution of final 
state mesons would also be of interest, as we could test the prediction of 
instanton-based isotropic production. 

In the second part of our work we turned our attention to the problem 
of nucleon spin. The experimental findings point to a large OZI violation in 
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the flavor singlet axial vector channel. Therefore we studied the OZI violation 
in axial-vector two-point functions and axial vector coupling of quark and 
nucleon. We found little reduction in the value of the flavor singlet coupling 
of quark, gg = 1, while some suppression was present in the coupling of the 
nucleon: g A = 0.77. However, the experimental value g A = (0.28 — 0.41) is 
significantly lower. Since the studied quark two and three point functions do 
not show large OZI violation we conclude that the structure of nucleon spin is 
given by another, possibly more complicated mechanism which does not show 
on the quark level, but is specific to nucleon. 

It would be interesting to study the connection to calculations of g A based 
on the matrix element of the topological charge density GG. One should also 
be able to get some insight from the lattice calculations of OZI violation in 
the axial vector channel. 

In the last part of present work we elaborated an algorithm for group 
integration and exemplified it on SU(N) fundamental and adjoint represen- 
tations. Although straight-forward, we hope the reader will find it a useful 
gadget, handy but inexpensive, with no reason not to have it in one's toolbox. 
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Chapter 6 



Appendix 



A Spectral Representation 

A.l Nucleon Two-Point Function 

Consider the euclidean correlation function 

njfOc) = (v a (o)f(x)), (A.i) 

where r) a (x) is a nucleon current and a is a Dirac index. We can write 

nf(x) = n x (x)(x ■ 7 ) Q/3 + n 2 (x)5 a ?. (a.2) 

The functions n 12 (a;) have spectral representations 

POO 

Ux(x) = / dspx(s)D'(^,x), (A.3) 
Jo 

poo 

n 2 (x) = / dsp 2 (s)D(^,x), (A.4) 
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where pi,2(s) are spectral functions and 



Tfl 

D(m,x) = Ki(mx), (A. 5) 

4n 2 x 



2 

m 



D'(m,x) = -—^K 2 (mx), (A.6) 



are the euclidean coordinate space propagator of a scalar particle with mass m 
and its derivative with respect to x. The contribution to the spectral function 
arising from a nucleon pole is 

PiO) = \X 2 N \5(s - m 2 N ), p 2 (s) = \\ 2 N \m N S(s - m%), (A. 7) 

where Xn is the coupling of the nucleon to the current, (0\rj\N(p)) = \nu(p), 
and 7Ti7v is the mass of the nucleon. It is often useful to consider the point-to- 
plane correlation function 



Kf(r) = j d 3 xUf(r,x). (A.8) 



The integral over the transverse plane insures that all intermediate states have 
zero three-momentum. The nucleon pole contribution to the point-to-plane 
correlation function is 



Kf{r) = 1(1 +74 )^|A^| 2 exp(-m w r). (A.9) 
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A. 2 Scalar Three-Point Functions 

Next we consider three-point functions. Before we get to three-point functions 
of spinor and vector currents we consider a simpler case in which the spin 
structure is absent. We study the three-point function of two scalar fields <j) 
and a scalar current j. We define 

U(x,y) = (<l>(0)j(y)<l>(x)). (A.10) 

The spectral representation of the three-point function is complicated and in 
the following we will concentrate on the contribution from the lowest pole in 
the two-point function of the field <fi. We define the coupling of this state to 
the field <fi and the current j as 

(O|0(O)|$(p)) = A, (A.ll) 
W)tf(0)|*(p)> = F(q 2 ), (A.12) 

where F(q 2 ) with q = p — p' is the scalar form factor. The pole contribution 
to the three-point function is 

n(x,y) = A 2 J d 4 z D(m,y + z)F(z)D(m,x - y - z), (A.13) 

where D(m,x) is the scalar propagator and F(z) is the Fourier transform of 
the form factor. In order to study the momentum space form factor directly it 
is convenient to integrate over the location of the endpoint in the transverse 
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plane and Fourier transform with respect to the midpoint 



d 3 x I dV" , (^)i(r/2,^(T,^) = _^_exp(-mr)F(g 2 ). (A.14) 



The correlation function directly provides the form factor for space-like mo- 
menta. Maiani and Testa showed that there is no simple procedure to obtain 
the time-like form factor from euclidean correlation functions [100]. 

Form factors are often parametrized in terms of monopole, dipole, or 
monopole-dipole functions 







F D {q 2 ) 


= F D (0)C 


F MD (q 2 ) 


= F MD (0)- 



my 



m 



2 ' 

l v 

% 2 

V 



Q 2 + m 2 v 



/ m 



Q 2 + m\ \Q 



rrtr 



(A.15) 
(A.16) 
(A.17) 



with Q 2 = —q 2 . For these parametrization the Fourier transform to euclidean 
coordinate space can be performed analytically. We find 

F M (x) = m 2 v D(x,m v ) (A. 18) 

F D {x) = m v (-^D'{x,m v )-D{x,m v f) (A.19) 

Fmd(x) = 2 1 2 2 < —a 2 (D(x,m 1 ) - D(x,m 2 )) 

m\ — mf [ — mf 

+ ^2 (^D\x,m 2 )+D{x,m 2 ))\ . (A.20) 

777/2 J 

We also consider three-point functions involving a vector current j^. The 
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matrix element is 

(H P ')\j,(o)Mp)) = ^nq 2 )- (A.21) 

The pole contribution to the vector current three-point function is 

n„(z, y) = A 2 J d 4 z D{m, y + z)z IM F'(z)D(m, x-y-z), (A.22) 

with F'(z) = dF(z)/dz and z^ = z^/\z\. For the parameterizations given 
above the derivative of the coordinate space form factor can be computed 
analytically. We get 

F' M (x) = m 2 v D'(x,m v ) (A.23) 
F' D {x) = -^D(x,m v ) (A.24) 

F'mdW = ?' m % ( 2 1 ^ (D'(x,m 1 )-D'(x,m 2 )) 
— m\ y m<2 — m{ 

+ |D(x,m 2 )}. (A.25) 
A. 3 Nucleon three-point functions 

Next we consider three-point functions of the nucleon involving vector and 
axial-vector currents. The vector three-point function is 

(U a VNN )f(x,y) = (v a (0)V^y)f(x)). (A.26) 

The axial- vector three-point function is defined analogously. The nucleon pole 
contribution involves the nucleon coupling to the current 77 and the nucleon 
matrix element of the vector and axial vector currents. The vector current 
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matrix element is 



(N(p>)\V*\N(p))=u(p>) 



-u(p), (A.27) 



where the form factors Fx 2 are related to the electric and magnetic form factors 



via 



G E (g 2 ) = Hq 2 ) + ^Hq 
G M {q 2 ) = F 1 (q 2 ) + F 2 (q 2 ). 



(A.28) 
(A.29) 



The axial-vector current matrix element is 



(N(p')\A a \N(p))=u(p' 



GA(q 2 h, + ^G P (q 2 )(p'-p) f 



l 5 -u(p) (A.30) 



where Ga,p are the axial and induced pseudo-scalar form factors. 

We are interested in extracting the vector and axial-vector coupling con- 
stants gy = Fx(0) and g& = Ga(0). In order to determine the vector coupling 
gv we study the three-point function involving the four-component of the vec- 
tor current in the euclidean time direction. For simplicity we take y = x/2. 
We find that 



(IW)f (x,x/2) = U VNN (r)5^ + n 2 ^(r)( 7 4) a/3 , 



(A.3i; 
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where x^ = (0,r). The two independent structures II are given by 



III 



VNN\ 



|A 12 



+ 



/ d y 

t - 2y 4 



}T + 2y 4 



2x 2 



2xi 

mD{x 1 )D'(x 2 )F 1 {y) 



mD'(x 1 )D(x 2 )F 1 (y) 



D {xi)D [x 2 — — 

X1X2 2m 



^vnn( t ) 



\Xn\ 2 I d 4 y 



xFi{y) + \y\ 



r z - Ay% + Ay 
Ax\X 2 



vf2 



■D'(x 1 )D'(x 2 )+m 2 D(x 1 )D(x 2 ) 



fib 1 , , , , 777/ 1 . 

—D'{x x )D{x 2 ) + —D'(x 2 )D( Xl ] 
X\ x 2 



2m 



where Fi t2 (y) are the Fourier transforms of the the Dirac form factors Fi^ 2 (q 2 ), 
and we have defined x\ = (y, r/2 + y 4 ) and x 2 = (—y, r/2 — y 4 ). 

In order to extract the axial- vector coupling we study three-point functions 
involving spatial components of the axial- vector current. We choose the three- 
component of the current and again take y = x/2 with x M = (0, r). We find 



(n^)f(x,x/2) = nWW(75) Q/3 +n 2 4iVJV (r)( 737 5)^+ni JVAr (r)( 737 475) Q/3 , 

(A.32) 

with 
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Xi x 2 



(A.33) 



G' P (y) \ 
2m J ' 

(A.34) 
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where Ga,p{u) are the Fourier transforms of the nucleon axial and induced 
pseudo-scalar form factors. 

These results are quite complicated. The situation simplifies if we consider 
three-point functions in which we integrate all points over their location in the 
transverse plane. The vector three-point function is 

J d 3 x J d 3 y(n VNN )f(T,x ] T/2,y) = ^(l + li rP\\ N \ 2 ex V (-m N T), (A.35) 

where gy = Fi(0) is the vector coupling. Note that the three-point function 
of the spatial components of the current vanishes when integrated over the 
transverse plane. The axial-vector three-point function is 

J d 3 x J d 3 y(U ANN )f(r,x; r/2, y) = y ((1 + 74)7375)°" \^n\ 2 exp(-m N r), 

(A.36) 

where g A = Ga{0) is the axial- vector coupling. In the case of the axial- vector 
current the three-point function of the time component of the current vanishes 
when integrated over the transverse plane. This is why we consider three-point 
function involving the spatial components of the axial current. 

A. 4 Phenomenology 

In Fig. 16.11 we show the nucleon pole contribution to the vector and axial- 
vector three-point functions. We have used the phenomenological values of 
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x [fm] 



G, andG„ 




Figure 6.1: Nucleon pole contribution to the vector (upper panel) and axial- 
vector (lower panel) nucleon three-point function. The solid line shows the 
complete results, the dashed line is the contribution from the F\ and Ga form 
factors only, and the dash-dotted line corresponds to a point-like nucleon. We 
have used a nucleon coupling constant A = 2.2 fm -3 as well as phenomenolog- 
ical values for the form factors and coupling constants. 
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Figure 6.2: Ratio of the phenomenological parameterizations of the axial- 
vector and vector three-point functions. We have added a short distance con- 
tinuum contribution to the nucleon pole terms. The curves are labeled as in 
the previous figure. Note that both the solid and the dashed line will approach 
qa = 1.25 as x — ► oo. Also note that the solid line is in very good agreement 
with the instanton calculation shown in Fig. 13.131 

the iso-vector coupling constants 



G E (0) = 1, G M (0) = 4.7, 

G A (0) = 1.25, G P (0) = ~t~9a- 



(A.37) 



We have parametrized Ge,v and Ga by dipole functions with my = 0.88 GeV 
and rriA = 1.1 GeV. The induced pseudoscalar form factor is parametrized as 
a pion propagator multiplied by a dipole form factor with dipole mass tua- 

We observe that at distances that are accessible in lattice or instanton 
simulations, x ~ (1 — 2) fm, the typical momentum transfer is not small and 
the correlation function is substantially reduced as compared to the result 
for a point-like nucleon. We also observe that the and Gp form factors 
make substantial contributions. Fig. 16.21 shows that the ratio of the axial- 
vector and vector correlation functions is nevertheless close to the value for a 
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point-like nucleon, qa! 9v — 1-25. We observe that the ratio of point-to-point 
correlation functions approaches this value from above. This is related to the 
fact that the axial radius of the nucleon is smaller than the vector radius. As 
a consequence, the point-to-point correlation function at finite separation r 
"sees" a larger fraction of the axial charge as compared to the vector charge. 
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B Instanton contribution to quark three-point 
functions 



In this appendix we provide the results for the traces that appear in the single 
instanton contribution to the quark three-point function. Our starting point 
is the expression 



/ peon \ 
\ r A/VQQ) 



NZNZ C A/V \ r ZMm ~r r mZM) 



(b.i; 



see equ. (I3.61J1 . Due to the Dirac structure of the non-zero mode part of the 
propagator, the NZNZ term is the same for both the vector and axial- vector 
correlation functions. It has 4 parts stemming from combinations of the two 
terms of non-zero mode propagator equ. 13.171 
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r NZNZ 2 2 



H(x, z, y)x a °y /3 
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(B.5) 



where 



-i 



H(x,z,y) = 
and the Dirac trace T± is defined as 
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where 2,3,... is short for /i 2 , /i 3 , . . .. The SzMS m term is easily seen to be 
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(B.8) 



with <^(x) = pj{j^^fx{x 2 + p 2 ) 3 / 2 ). The S m SzM term is obtained similarly 



r mZM 
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(B.9) 
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For all the above formulas, we need to add instanton and anti-instanton contri- 
bution and integrate over the position of the instanton, which was suppressed 
above. As usual, for an instanton at position zj we have to shift x, y, z accord- 
ing to x — > (x — zj), etc. 

B.l Computation of path exponent 

The correlation functions of operators at different points are made gauge in- 
variant by introducing the connector, the Wilson line. In this appendix we 
compute the Wilson line in instanton background. 

In order to perform the computation, we need to evaluate the path expo- 
nent 



For a general gauge configuration that is not an easy task as the path exponent 
is defined as follows. Parametrize the path by s G [0, 1]. The integral in the 
exponent is then a limit of a sum of matrices, which might not commute. 
Expand the exponential in power series and reorder the matrices in each term 
in decreasing value of parameter s. The sum obtained this way is the well- 
known Wilson line. 

For the path exponent along a straight line in the field of an instanton 
the computation is actually much simpler. Let us parametrize the line by 
x^ = x^ + s(x2 — Xi)^, so that the end points Xi and 22 correspond to s = 




(B.10) 
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and s = 1 respectively. The integral then reads: 



pX 
J X\ 
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L 
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p 2 x v (s) dx^ , 

ds 

x 2 (s)(x 2 (s) + p 2 ) ds 

P 2 [xi{x 2 - xiY + s(x 2 - xiY(x 2 - XiY] 
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AYT a dx» 



n- a t' 



.a 



1 




M 



ds, 



with r a Pauli matrices. 

The second term in the numerator vanishes due to the anti-symmetry of 
t'Hooft's tensor rf = — 77" . Then the integrand is in fact a constant matrix 
multiplied by a function of s at any point on the line. As a result, no path- 
ordering is necessary since the matrices at any point commute. Therefore we 
can integrate first and then expand the exponential to get: 




Xl 




O2 - 



x 



{fli(xi,x 2 ) [tan -1 (x 2 ■ (x 2 - X\)Q,\{xi, x 2 )) 
-tan" 1 (xi ■ (x 2 - Xi)Qi(xx,X2))] 

Q 2 (xi, x 2 , p) [tan" 1 (x 2 ■ (x 2 - x 1 )Q 2 (x 1 ,x 2 , p)) 

- tan" 1 (xi • (x 2 - Xi)Q 2 (xi, x 2 , p))] } (B.ll) 



where 
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The expansion of the exponent of Pauli matrices is then easy: 



Pexp 



(i JJ A a ^T a dx»j = cos[|iir ± (xi,z 2 )|] 

+i sin [\K±{x u x 2 ) |] n x {x x ,x 2 )ff a x\{x 2 - x x Yr a (B.12) 



where the modulus of K a (xi,x 2 ) is 



\K(xi,x 2 )\ = VK a K a 

= {tan^ 1 (x 2 ■ (x 2 - x 2 )) 

— tan -1 (xi • (x 2 — Xi)Qi(xi,x 2 )) 



Q 2 (x 1 ,x 2 ,p) 



(tan 1 (x 2 ■ (x 2 - Xi)Q 2 (xi,x 2 ,p)) 



- tan -1 (xi ■ (x 2 - Xi)Q 2 (xi,x 2 , p))) } (B.13) 



In the above formulas we took the instanton to be centered at 0. The general- 
ization to instanton center at z is straight-forwardly obtained by the replace- 
ment Xi — > Xi — z and x 2 — > x 2 — z. 

For example, for the case of an infinite line in 4-direction we get: 



+oo 



Pexp ( % I A a /r a dx 4 



cos 



% sin 



nil 



nil 



\A ■ 



t (B.14) 



in agreement with formula (206) of [24]. 
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C Euclidean matrices, conventions, 

We use the following conventions: 



gT = 6^ = diag(l,l,l,l) = 9^ 



,1234 



£1234 



The contractions of e therefore are 











£pi>af3^^ 


K 















(c.i; 



and e pp , Q/3 e^> = 2(<^ - 5"J%), e pp , a(3 e^'^ = 65% and e Pf/aP ^ = 24 

Euclidean gamma matrices are defined as {Y>Y} = 2g^ v . We take Y 
^ 1 ^ 2 ^ 3 ^/ 4 a nd <T Miy = \[Yi 7"]. The following relations apply: 7^ = 7^, 75a 
_l e ^ua/3 . Q ne can eas jiy compute the traces: 



pU 



Tr[YY] = 4g 

Tr[YYYY] 

Tr[YYYYi 5 } 

Tr\YYYY' Y' Y\ = 4 
Tr [YYYY'Y'Yi 5 ] = 4 



A(g^g u P - g^gf" 7 + g" p g au ) = 4S pavp 

^povp 

gP<J gpp'a'v _ gppgap'a'v _|_ /jp' gapa'v 
_gH<j' gvpp'v _|_ gPVgapp'a' 



1 ,,/ _/_/. 



^Mjpv Qn>' p'cr'v _|_ gpapv p a 



(C.2) 



Some useful decomposition: (for transparency we will short our notation from 
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7 V7 3 7 4 7 5 7 6 7 T = 7 T (/ M + ^)^[^2,3,4,5,6] (C.3) 

where the Dirac trace was defined in IB. 71 The four vector of SU(2) C 
matrices is o± = (a, =pi) while the color - rotated matrices are given as cr±(R) = 
(R ab a b , =p). Here R ab is the adjoint representation rotation matrix, defined as 
R ab a b = u a au^ wit h u e SU{2). The following relations link the SU{2) 
generators to t'Hooft r\ tensor: 

o%<?_ = + i v a ^R ab cT b 

oto" + = g" u + if] a ^R ab a b (C.4) 

The definition and useful properties of rj can be found e.g. in [24]. For our 
computations suffices to list the most often used formulas: 

n%-rfp\ = 5 W $v\ ~ ^xKp T £pup\ (C.5) 
£ abC ?faf3 ft h pu~rfp\ = ^p'faP ^fvX ~ ^pXffa/3 ^Ip +^A^f^ p ffa/3 ~ ^p'T^X^aff 

where -f) is fj (77) for upper (lower) sign. 

After tracing the color indices, all the rotation matrices cancel due to 
orthogonality R ab R cb = 5 ac and unimodularity det R — 1. The following traces 
are therefore also valid for color-rotated matrices. 

Tr[a^a v _} = Tr[a^a u + ] = 26^ (C.6) 
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Tr[a^a p ± a^] = 2{S" vpa T e^ pa ) (C.7) 

As one can easily see, the traces of SU (2) matrices Tr[a± ■ ■ ■ ] can be obtained 
from the traces of corresponding 7 matrices projected to the particular 2x2 
submatrix as follows: Tr c [a± ■ ■ ■ aV[ = Ttd[{^ ■ ■ ■ 7^)7 T ] 

We also provide some useful decompositions for the multiplication of color 
matrices: 

°T a ± a T a ± = 5 p<T 5^-^ a ' P V^+^ a L a ^S pu +^ a '^5 pa ^ b ' p(7 ^ c ' pu e bca (C. 



*v „Oi „(3 „p (7 



1^6 pv 5 a ^ 5 pcJ — s fil/ ff a,al3 ^f ,pc7 

_ ^ c ^ v §P a _|_ fj. c 'P a S a/3 — ?f a > a P ff b ' pa e abc I 

_ e ecf ^ e >M^ ^c,a/3 ^pa _|_ ^c,pa _ ^_a,a/3 ^a! ,pa e aa'cj| 



(C.9) 
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